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Abstract 

This collection of studies on new physics at the LHC constitutes the report of 
the supersymmetry working group at the Workshop 'Physics at TeV Colliders', 
Les Houches, France, 2007. They cover the wide spectrum of phenomenology 
in the LHC era, from alternative models and signatures to the extraction of 
relevant observables, the study of the MSSM parameter space and finally to 
the interplay of LHC observations with additional data expected on a similar 
time scale. The special feature of this collection is that while not each of 
the studies is explicitely performed together by theoretical and experimental 
LHC physicists, all of them were inspired by and discussed in this particular 
environment. 
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Aim and structure of this collection 1 

With the first LHC data around the corner, the great common goal of theoretical and experimental high- 
energy physics appears to be in close reach: over the coming years, we will have to try to understand 
the origins of electroweak symmetry breaking and the role of the TeV scale from the data rolling in. 
This effort can only be successful if theorists and experimentalists work in close collaboration, following 
the well-known spirit of the Les Houches workshops. This collaboration of course starts with the proper 
understanding of QCD, the theory which describes any kind of particle production at the LHC, but which 
also describes the main backgrounds which Higgs and new-physics searches have to battle. However, 
due to the complexity of LHC data on the one hand, and of the new-physics models at the TeV scale on 
the other hand, the interaction between theorists and experimentalists needs to go much further. Realis- 
tically, we expect that any new-physics search at the LHC will require theorists to formulate viable and 
predictive hypotheses which are implemented in state-of-the-art simulation and extraction tools. Such 
models can guide experimental searches towards the of course yet unknown ultraviolet completion of the 
Standard Model, even if at least all but one known models for new physics at the LHC will be soon ruled 
out. 

Over recent years, high-energy theorists have hugely expanded the number of viable ultraviolet 
completions of the Standard Model. The main guiding principle of all of these models is still electroweak 
symmetry breaking. There are essentially two paths we can follow to explain the weak-scale masses of 
gauge bosons and of (third-generation) fermions: first, we can assume the minimal Higgs sector of 
our Standard Model to hold, which leads to the hierarchy problem. Without solving this problem, the 
Standard Model appears to be incomplete as a fundamental theory valid to all mass scales up to the 
Planck scale. Such ultraviolet completions are particularly attractive if they allow us to incorporate 
dark-matter candidates or unification scenarios. Supersymmetry is in particular in the experimental 
community definitely the most carefully studied completion, but extra dimensions or little-Higgs models 
are alternatives worth studying. The alternative path to describe electroweak symmetry breaking are 
strongly interacting models, which avoid predicting a fundamental Higgs boson. Such models have 
recently become more viable, if combined for example with extra dimensions. 

As indicated by its title, this working group focuses on supersymmetry as one example for an 
ultraviolet completion of the Standard Model with the usual Higgs mechanism. However, simply writing 
down one supersymmetric version of the Standard Model does not suffice in view of the almost infinite 
number of LHC analyses which would be possible to study such models. As a matter of fact, signatures 
which until recently were thought to be typical for supersymmetry, namely jets plus missing energy plus 
maybe like-sign dileptons are by now mainstream signals for extra dimensions, little-Higgs models, 
or even strong interactions. Therefore, this collection of projects should first be considered as studies 
of models which lead to typical supersymmetry signatures, mostly beyond the naive inclusive 'missing 
energy plus jets' analysis. The obvious question is how with the LHC running we would go about to 
understand what the underlying theory of such signatures could be. Secondly, we include studies on 
version of a supersymmetric Standard Model which deviate from the simple MSSM. Which means that 
even supersymmetry as an underlying principle does not have to look exactly like we naive think it 
should look. It is a healthy development that theoretical physics has moved beyond its focus on the 
minimal supersymmetric ultraviolet extension of the Standard Model, while at the same time, not all 
alternatives in and beyond supersymmetry need to be studied including full detector simulations. 

Independent of our 'preferred' models, theorists need to carefully formulate viable TeV-scale 
models, including variations of key models which allow us to test predictions for example of the MSSM. 
These alternatives can for example be driven by the similarity of signatures or by the aim to test certain 
underlying theory structures. In supersymmetry such models obviously involve Dirac gauginos (altering 
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the Majorana nature for example of the gluino or the dark-matter agent). Naively matching the two 
degrees of freedom of an on-shell gluon with those of a gluino only allows us to write down a Majorana 
gluino. In extended supersymmetric models this requirement does not need to exist. On the other hand, 
even assuming minimal supersymmetry an extended Higgs sector like in the NMSSM, should be tested 
carefully. Moreover, supersymmetric spectra with a hierarchy between gauginos and scalars maintain 
many advantages of supersymmetry, like unification and a dark-matter agent, while avoiding flavor con- 
straints at the expense of introducing fine tuning. Due to the organization of the complete document, in 
this chapter of the Les Houches proceedings we limit ourselves to variations of the MSSM, deferring 
models like extra dimensions or little-Higgs to another collection of articles in the same volume. 

In a second step, theorists and experimentalists need to develop strategies to extract information on 
TeV-scale physics from LHC data. In the standard MSSM scenarios, studying the kinematics of cascade 
decays has been shown to be a spectacularly successful, even in the presence of missing energy from a 
dark-matter agent escaping the detector unseen. There are, however, many more or less experimentally 
complex ways to use LHC observables to extract information on the masses of new states from LHC data. 
The information can come from the general underlying mass scale is event samples including physics 
beyond the Standard Model, or from any combination with cascade information. Of course, such studies 
are not limited to supersymmetry, but they can in principle be used for any new-physics signatures with 
decays from strongly interacting new-physics states down to a weakly interacting dark-matter particle. 

Whatever we are looking for at the LHC, technically correct simulations of new-physics events 
are crucial, if we ever want to extract the fundamental parameters from their comparison with data. 
There is no good reason to try to extract new physics from 21st-century LHC data using 20th-century 
Monte-Carlos and methods. The past years have seen impressive progress in incorporating new-physics 
signals in modern Monte-Carlo tools, including for example the proper simulation of many-particle 
final states beyond a naive narrow-width approximation. This particular problem is being studied in 
the supersymmetric framework, but it is at least as relevant for new-physics models which predict more 
degenerate mass spectra, like for example generic universal extra dimensions. 

Obviously, LHC data on TeV-scale physics will not come into a data-free world. There is a 
wealth of information we have already collected on such physics models, and during the LHC era we 
expect much more of it. The long list of current and future complementary data includes electroweak 
precision data, the muon's anomalous magnetic moment, precision-flavor physics, dark-matter mea- 
surements linked with big-bang nucleosynthesis, and most importantly at some stage the high-precision 
data from a future ILC. In particular, when it comes to measuring as many model parameters of the 
TeV-scale Lagrangian as possible, the proper combination of all these pieces of information is crucial to 
our understanding of the ultraviolet extension of the Standard Model. Only once we can claim a solid 
understanding of the TeV scale we should attempt to extrapolate our physics picture to very large mass 
scales, to finally determine if our underlying theory can really be a fundamental theory of Nature. After 
all, the LHC is not going to be the last, but the first major experiment allowing us to carefully study the 
TeV scale and determine the fundamental parameters of physics beyond the weak scale over the coming 
years! 

The last section in this collection follows a great tradition of the Les Houches workshops: the 
successful definition and implementation of interfaces between computer tools used by the theoretical 
and experimental high-energy community. 

At this stage, the conveners of the SUSY session would like to express their gratitude of course to 
the organizers of this inspiring and enjoyable workshop. Moreover, we would like to thank all the young 
collaborators in Les Houches and elsewhere, who have made possible the impressive studies presented 
in this collection. 
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1. A model for dirac and majorana gaugino masses 1 
1.1 Introduction 

Massive fermions can appear either as Majorana or Dirac. Because the latter allow charged states, they 
are easier to detect . And, in fact, all identified fermionic masses are of Dirac type. The nature of neutri- 
nos masses remains unknown, and unveiling it is the main challenge for double beta decays experiments. 
It is then legitimate to ask about the form of the masses of new fermions that could be detected by LHC. 

In the minimal extension (MSSM) masses of gauginos are of Majorana type. Obtaining Dirac 
ones requires pairing up with new fermions that should then arise as components of extra chiral fields in 
the adjoint representation. The easiest way to incorporate these new states is to make the gauge fields as 
parts of A = 2 multiplets. Such a scheme is present in an extra-dimensional picture where the N = 2 
fields appear as bulk states (in the absence of a projection) while the chiral matter appear as localized 
states in N = 1 representations. In such a set-up Dirac masses have been shown to appear naturally 
in the presence of an anomalous U (1) as a result of new operators that mix the MSSM fields with the 
anomalous U(l) [1, 2, 3, 4, 5, 6]. It was later shown that such operators can be generated by loop 
effects if the supersymmetry breaking sector is in an N = 2 representation [7, 8, 9], and that the quartic 
tree-level Higgs potential is in fact also modified [9]. Such models where studied as primarely based on 
the only presence of D-term breaking suffer from two issues. First, in the minimal set-up the adjoint 
scalar have tachyonic masses[2, 9]. Second, typical supersymmetry breaking model would lead to both 
D and F-terms. The latter will fix the first issue, but also turn on new sources for sof terms, in particular 
Majorana masses for the gauginos. 

Here, we will take a different, more phenomenological approach. We will provide with the corre- 
sponding Lagrangian containing an N = 2 extended gauge sector, and we will not address the origin of 
the soft masses. 



1.2 Primer 

We will start by fixing our conventions for the spinorial notation as well as for the MSSM, and then we 
will proceed to extend the gauge sector. 



Spinors notation 

A Dirac fermion \I> d has four components which can be assembled into two-component spinors : 

MS) 

where ifj a = e Q/3 V ( g an d i>a. = (V'a)*- Here e af3 is the completely antisymmetric tensor and e 12 = 1. We 
adopt for the 7 matrices the following representation: 

v = (° ;'')■ t" = *VtV = (J 

where 7r are the Pauli matrices and the notation cj m = (1, ~a) and a^ 1 = (1, —It) is used. A Dirac mass 
term here takes the form: 

~*d^d = xV> + Y>X, (3) 

with the spinor products xV' an d xY' defined as = X a "4 ) a and xY' = XdV'"' where, again, we use the 
notation \ a = e Xf3- A Majorana fermion can also be written as: 
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with the Majorama mass term: 

*m*m = XX + XX, (5) 
Below we will always use left handed fermions. The Dirac fermion representing a lepton is: 

*'-(£!) 

where %^ ) is the left-handed lepton field and x^ l+ ^ is the charge conjugated of the left-handed anti- 
lepton field: = £ a ^(Xa ^)*- Then, its Dirac mass will be written molx^ ^X^ l+ ^ + X^ + ^X^ 

Super symmetry 

The generic supersymmetric Lagrangian density for a gauge theory discussed here take the form 2 . 

= £ gauge ~i~ ^--chiral ~\~ ^--minimal coupling- (7) 

Here, the gauge kinetic Lagrangian is given by: 

Cgauge = ~F^F^ a + iT<T^D u \) a + ^D a D a (8) 

where F uu is the gauge boson field strength, A and D a are the associated gaugino and and auxiliary field, 
respectively. is the gauge covariant derivative. Here, the index a is a gauge symmetry group index 
corresponding to the generator T a . 

The chiral Lagrangian is written as: 

fdW 1 d 2 W \ 
Cchiral = D^D^i + ixfDnXi + F*Fi + I — Fi - --^-^- XiXj + h.c.) . (9) 

Here the chiral fermion Xi , the boson fa, and the auxiliary field Fi belong to the same gauge group 
representation and form an N = 1 multiplet . The index i labels the different chiral multiplets. The 
superpotential W is an holomorphic function of the fields 0j. 

Finally, the last piece in the supersymmetric Lagrangian density is: 



Zmimmal coupling = -g{4>*T a ^)D a - V2g feT^A" + A°(XiT°&) 



(10) 



where g is the gauge coupling contant. The last two terms in ( 10) will be important to us as the scalar 
field takes a vacuum expectation value (the Higgs multiplets), producing bilinears in the fermions, thus 
mass terms. 



1.3 MSSM 

The field content of the MSSM is presented in table 1. Note that all the chiral fermions are left-handed, 
the charge conjugation label c allows to use the appropriate antiparticles. At the renormalisable level, the 
MSSM has the superpotential: 

W = y i iu c i Q j ■ H u - y'JdlQj ■ H d - y^Lj ■ H d + fiH u ■ H d . (11) 

The indices i,j are family indices and runs from 1 to 3. The 3 x 3 y matrices are the Yukawa couplings 
and the parameter \i is a Dirac mass for the Higgsinos. The "■" denotes the SU(2) invariant couplings, 
for example: Q ■ H u = ulH® — cLlH+. 



2 In this subsection and the next we follow closely the presentation of Refs [10]. 
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Names 




Spin 


Spin 1/2 


Spin 1 


SU(3), SU{2), U(1) Y 


left-handed 


Q 


(uL,d L ) 


(u L ,d L ) 




3, 2, 1/3 


quarks 


u c 


u\ 


u l 




3, 1, -4/3 


(x3 families) 


d c 


d°L 






3, 1, 2/3 


leptons 


L 




{v eL ,e L ) 




1,2,-1 


( x 3 families) 


c 

e 


e L 


r 

e L 




1, 1,2 


Higgs 


H u 








1,2,1 




H d 








1,2, -1 


gluons 


g 




9 


9 


8,1,0 


W 


W 




W ± ,W° 


W ± ,W° 


1,3,0 


B 


B 




B 


B 


1,1,0 



Table 1 : Chiral and gauge multiplet fields in the MSSM 



Soft breaking 

The breaking of supersymmetry is parametrized by a set of terms, labelled soft as they preserve the 
absence of quadratic divergences. The possible soft breaking terms in the MSSM are quite limited, there 
are gaugino masses for each gauge group, squarks mass terms, sleptons mass terms, Higgs mass terms 
and triple scalar couplings. We are primarily interested the gaugino masses given by: 

-^{M 3 g a g a + M 2 W a W a + M l BB + h.c.) (12) 

Note that these are Majorama masses. Two of these terms for W 1 ' 2 combine as a Dirac mass for W ± . 

Neutralino masses 

The neutral fermions of interest are the higgsinos H® and and the gauginos, bino B and wino W°. 
The mass terms for these fields in the MSSM have three origins: 

• The soft breaking terms for the bino and wino, 

--(M 2 W°W° + M 1 BB + h.c). (13) 

• The two last terms in equation (10) generate a mixing between (B,W°) and (H^,H^). These 
mass terms are parametrized by the vev of the Higgs scalars < H® >= v u and < >= v d ■ 
Using tan/3 = v u /v d , one can then express them in terms of (3, the weak gauge bosons masses 
mw and mz and the weak mixing angle 9w as 



-m z 



cos 6 l w {cos H^W° - sm PH°W°) + sin 6 w (sin (3 H°B - cos B) + h.c. (14) 
• The n term in the superpotential W contributes to the higgsinos masses , 

tiH°,H° d + h.c. (15) 

Chargino masses 

Here we consider the charged higgsinos H+ and and the charged gauginos W + and W~. In the 
MSSM the origin of the chargino mass terms is completely analogous to those presented in subsection 
1.3, they take here the following form: 

-M 2 W + W~ (16) 
-V2m w s\n(3H+W- - V2m w cos (3H~W + + h.c. (17) 
-liH+,Hj + h.c. (18) 
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1.4 Extended susy gauge sector 

We now consider the scenario where the gauge sector arise in multiplets of N = 2 supersymmetry 
while matter states are in N = 1 SUSY representations. Moreover, the Higgs multiplets H u and are 
assumed to form an N = 2 hypermultiplet. 

The field content for the gauge sector is described in table 2. Note that for each N = 1 gauge 
multiplet present in the MSSM one need to introduce one extra scalar and fermionic fields. The latter are 
differentiated by a symbol ' (see table 2). 



Names 


SpinO 


Spin 1/2 


Spin 1 


SU(3), SU(2), U(1) Y 


gluons 

W 


y± yO 


9,~9' 

W ± ,W°,W' ± ,W'° 


9 

W ± ,W° 


8,1,0 
1,3,0 


B 




B, B' 


B 


1,1,0 



Table 2: N = 2 gauge supermultiplets fields 



In addition to the Majorana masses, Dirac ones can now be written. We will extended the MSSM 
soft terms to include 



1 



(M^g' a g' a + M' 2 W ,a W ,a + M[B'B') - {M^~g a ~g ,a + M?W a W ,a + M^BB') + h.c. (19) 



rD-fira-rirfa 



rD - 



where Mj are Majorana and are Dirac masses. 

The N = 2 supersymmetry in the gauge sector introduces new couplings analogous to the two last 
terms in equation (10). These lead to new bilinear mixing terms between gauginos and higgsinos when 
the Higgs scalars H® and H® acquire vevs. 

• Neutralinos: 



-mz 



sin 9 W (sin/J^B' + cos f3H^B') - cos e w (cos f3H^W /l) + sin fiH^W'") + h.c 



• Charginos: 

- V2m w cos f3H+W'- + V2m w sin (3HJW >+ + h.c. 
1.5 Fermionic mass matrix 



(20) 
(21) 



We now put all the previous terms together and describe the resulting mass matrices for both neutral and 
charged gauginos and higgsinos when both Majorana and Dirac term are present. 



Neutralinos 

The neutralino mass terms are presented in equations (13 )-(15), (19) and (20). The neutralino mass 
matrix, in the (B',B, W'°, W°,H^,H^) basis is: 






M N{ 



eut 



( M> 


Aff 




Mi 
















-mzswcp 


\ m z s w cp 


mzswsp 


cos 9w, sw = 


sm#jy, C/3 



M' 2 



-mzcwcp 






mzswsp 


mzswcp 


\ 




Mf 


-mzswcp 


mzswsp 




-mzcwsp 


-mzcwC/3 




M 2 


mzcwcp 


-mzcy/sp 




mzc w cp 





-V 


/ 


-mzcw-sp 








(22) 



Clearly this 6x6 matrix will provide in general very long expressions for the neutralino mass 
eigenstates. A simple case is when the mz dependent terms in (22) are relatively small compared to 
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the other entries, and can be treated as perturbations. Moreover, the gaugino Majorama masses are 
symmetric in the primed and unprimed fermions: M{ = M\ and M' 2 = M 2 . In this case the higgsino 
mass eigenstates are given (approximately) by the combinations 

n° 8 ~±{H° + H°) , JBfl~-L(fl2-fl2) (23) 

both having mass squared ft 2 . The neutral gaugino mass eigenstates are given, to leading order in 

mz/Mi, by 

B S = ±=(B + Bi) , B A = ±={B-&) (24) 

W° s = -^(W° + W'°) , W° A = -^(W° - W'°) (25) 
V2 V2 



with masses 



m §s ~M 1 + Mf ) , m SA ~ Mi - (26) 
m W o ^ M 2 + M 2 D m W o ~ M 2 - M 2 D . (27) 

One can express the ratio between Dirac and Majorama masses by the angle 6 D defined by tan 8 D = 
M D /M. Alternatively, the angle 6 D is measured by 

2 2 
m~ — m~ 

sin20g = -f2 §± (28) 

m- + m~ 

and 

2 2 
n m VKg ~~ m W° A 

sm29& = -^ p. (29) 

Charginos 

The mass terms for the charginos can be expressed in the form 

-\{{v~f 'M Ch v + + (v+fM^iT + h.c) (30) 

where we addopted the basis v + = (W' + , W + ,H+), v~ = (W'~, W~, H d ). Collecting all the terms 
presented in equations (16)-( 18), (19) and (21) leads to the chargino mass matrix : 

/ M' 2 M 2 D V2m w cos/3 \ 

M Ch =\ M 2 D M 2 y/2m w s\nf3 . (31) 

\ — y/2mw sin (3 y/2mw cos (3 fi j 

This nonsymmetric mass array is diagonalized by separate unitary transformations in the basis v + and 
v~ , M^ 9 = U^MchV, where the matrices U and V are unitary 

For the simple case considered in the previous subsection, where the mz dependent terms can be 
treated as perturbations and M 2 = M 2 , the higgsino mass eigenstates will be given approximately by 
and HT, both having mass squared fi 2 . The charged gaugino mass eigenstates will be given, to 
leading order in mz/Mi, by the combinations 

Wg ~ -^=(W + + W' + ) , W\ ~ -^=(W + -W' + ) (32) 
V 2 v 2 
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W s ~ -WW" + W'~) , WX ^ —(W- - W'~) (33) 



with squared masses 



~ (M 2 + M 2 D ) 2 m 2 ~ (M 2 - M 2 D ) 2 (34) 



~ (M 2 + Mf) 2 ~ (M 2 - M 2 ^) 2 (35) 



Note that, in this limit, the winos have approximate degenerate masses: m^ rs ~ m ^ 7 anc ^ 



Gluinos 

Since gluinos 5 and g' are in color octet representation, they cannot mix with any other fermion, the only 
possible gluino masses are the soft ones presented in (12) and (19). In the basis (<?', g) the gluino mass 
matrix is simply 

( ML MP \ 

m ^ = uj> m 3 ) ■ (36) 

We will illustrate two limits. The first one is when gaugino Majorama masses are symmetric in the 
primed and unprimed fermions: M3 = M3. The analysis of the gluino mass matrix in this case follows 
closely the discussion after (23). The gluino mass eigenstates are 

1 1 

9s = —^(9 + 9') , 9A = —^(9 - 9') (37) 

with masses 

m- gs = M 3 + MP, m u = M 3 - MP, (38) 

and the ratio between Dirac and Majorama masses is parameterized by the angle t&nOP = MP /M3, 
alternatively by 

2 2 

sin2#f = — p §±. (39) 

m 9S + m 9A 

The second limit is when one of the Majorama masses, say M' z , is very small compared to the 
other entries of the mass matrix (36). In this limit, the gluino mass eigenstates are given approximately 
by 

gi ~ cos a g — sin a g' , £f 2 — sin a g + cos a g (40) 

where 

tan a = -- cot OP (l + y/l + 4 tan 2 0f ) (41) 
and, as before, tan 6^ = MP/M3. The gluino masses are 



m 5l ~ T (l + ^l + 4tan 2 ^) , m §2 ~ — (l - y/l + 4 tan 2 OP) (42) 
and the ratio between Dirac and Majorana masses is parameterized 

sin 2 ef = 2 Jl^gpJ (43) 
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1.6 Interactions 

We turn now to the interactions between the extended supersymmetric sector and the MSSM fields. First, 
we remind the interactions between the MSSM gauginos, the higgsinos and the scalar Higgs: 



9 

V2 



w 



+ HWH d W l ) ~ ^= {H* U H U B - H* d H d B) 



(44) 



where g and g' are the SU{2) and U (l)y coupling constants. 

Due to the fact that the two Higgs H u and H d form a N = 2 hypermultiplet, their interactions 



with the new fermions W' and B' are 

H u ■ (a l H d )W' 1 + H d ■ (a l H u )W'^ -j=(n d - H U B' - H u ■ H d &) . (45) 
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One can straightforwardly verify that these interactions lead upon electroweak breaking to the gaugino- 
higgsino mixing present in the neutralino and chargino mass matrices . 



1.7 New scalars 

The N = 2 vector multiplets include, in addition to the new fermions dicussed above, scalar fields in the 
adjoint representation. We denote these states as T, g , T,\y and They couple to the Higgs chiral fields 
(now in an N = 2 hypermultiplet ) in the superpotential, and through their F-term modify the tree-level 
Higgs scalars quartic terms in the potential by the new terms: 

-y E \ H - ■ al M 2 ~Y lHu ' Hdl2 ■ (46) 

i 

These scalar should not remains massless, but get soft terms as: 

with mj s > mj A . If the masses in (47) are big compared to the Higgs mass, these fields can be integrated 
out and in the low energy theory the scalar potential is the one in the MSSM plus the contributions coming 
from the quartic terms (46) [9]. Note that the latter contribution disappears, if instead, the integration out 
is supersymmetric (due to a large supersymmetric mass). 
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2. NMSSM in disguise: discovering singlino dark matter with soft leptons 1 
2.1 Introduction 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) provides an elegant solution to the fi 
problem of the MSSM by the addition of a gauge singlet superfield S [1 1, 12, 13, 14]. The superpotential 
of the Higgs sector then has the form XS(Hd ■ H u ) + ^kS 3 . When S acquires a vacuum expectation 
value, this creates an effective fi term, [i = X(S), which is automatically of the right size, i.e. of the order 
of the electroweak scale. 

The addition of the singlet field leads to a larger particle spectrum than in the MSSM: in addi- 
tion to the MSSM fields, the NMSSM contains two extra neutral (singlet) Higgs fields - one scalar and 
one pseudoscalar - as well as an extra neutralino, the singlino. Owing to these extra states, the phe- 
nomenology of the NMSSM can be significantly different from the MSSM; see chapter 4 of [15] for 
a recent review and references. In particular, the usual LEP limits do not apply to singlet and singlino 
states. Moreover, the singlino can be the lightest supersymmetric particle (LSP) and a cold dark matter 
candidate. 

In this contribution, we investigate the LHC signature of a SPS la-like scenario, supplemented by a 
singlino LSP. In such a setup, gluinos and squarks have the 'conventional' SUSY cascade decays into the 
bino-like neutralino, x\ ~ which then decays into the singlino LSP, Xi ~ S, plus a pair of opposite 
sign same-flavour (OSSF) leptons. (The x% decay proceeds dominantly through an off-shell slepton.) A 
dark matter relic density of Qh 2 ~ 0.1 is obtained if the Xi and/or annihilate through pseudoscalar 
exchange in the s-channel. 

One peculiar feature of this scenario is that the mass difference between X\ and x\ * s always small; 
it reaches at most ~ 12 GeV, and is often much smaller. The leptons originating from the x% ~^ Xil + l~ 
decay hence tend to be soft. In the standard SUSY analysis, requiring pr{l^) > 20 GeV, there is a 
risk of missing these leptons and wrongly concluding to have found the MSSM instead of the NMSSM, 
with x\ as the LSP and dark matter candidate (discovery of the additional Higgs states will also be very 
difficult at the LHC in this scenario). The aim of this contribution is to show the feasibility of detecting 
the x\ Xi^~ decay and measuring the singlino-bino mass difference by looking for soft di-leptons. 

We use the NMHDECAY [16, 17] program to compute the NMSSM mass spectrum and Higgs 
branching ratios, and to evaluate the LEP bounds; SPHENO [18] is used to calculate the sparticle 
branching ratios, and MICROMEGAS [19, 20] for the relic density. The SUSY-breaking parameters of 
our scenario are listed in Table 1. The main difference to the original SPS la [21] is that we choose 
Mi = O.5M2 = 120 GeV, leading to a X2 mass of ~ 115 GeV, in order to evade LEP bounds when 
adding the singlino and singlet Higgses. To obtain a singlino LSP, we choose A ~ 10~ 2 and k ~ 0.1 A. 
This way Xi ~ 99% S, and m^o hardly varies with A and k (~ 0.1 GeV). In addition, the trilinear Higgs 
couplings A\ and A K are chosen such that m^o + m^o ~ m^, for at least one combination of i, j = 1, 2, 

Xi Xj 

to achieve 0.094 < VLh? < 0.135 [22]. We thus obtain a set of NMSSM parameter points with varying 

Am = m^o — m^o . 

The four points used for this study are summarised in Table 2. Points A-D have Am = 9.7, 3.0, 
1.5 and 0.9 GeV, respectively. The SM-like second neutral scalar Higgs, S2, has a mass of 115 GeV for 
all these points, consistent with the LEP limit. On the other hand, the lightest neutral scalar Si and the 

l S. Kraml and A.R. Raklev 



Mi M 2 M 3 /ieff M iis M ti Mq i M &1 

120 240 720 360 195 136 133 544 526 524 496 420 521 



Table 1: Input parameters in [GeV] for our SPS la-like scenario. The NMSSM-specific parameters are given in Table 2. 
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Point 


A[1(T 2 ] 


K[W- :i ] 


A x 


A K 




m Al 


m A2 


m Sl 
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A 


1.49 


2.19 


-37.4 


-49.0 


105.4 


88 


239 


89 


0.101 


7 x 10" 


ii 


B 


1.12 


1.75 


-42.4 


-33.6 


112.1 


75 


226 


100 


0.094 


9 x 10" 


-13 


C 


1.20 


1.90 


-39.2 


-53.1 


113.8 


95 


256 


97 


0.094 


1 x 10- 


-13 


D 


1.47 


2.34 


-39.2 


-68.9 


114.5 


109 


259 


92 


0.112 


4 x lO- 


-14 



Table 2: NMSSM benchmark points used in this study. Masses and other dimensionful quantities are in [GeV]. 




P T [GeV] 

Fig. 1: pr distributions for leptons from the decay X2 ~* Xil + l~ m benchmark points A-D. All distributions are normalised 
to unity over the whole momentum range. 



lighter pseudoscalar A\ are mostly singlet states, and can hence be lighter than 1 14 GeV. Concerning the 
neutralino annihilation, for Point A the dominant channel is X2X2 — * contributing 88% to (av). For 
Point B, XiXi> X1X2 an d X2X2 annihilation to bb contribute 10%, 15%, and 50%, respectively. Point C 
has again dominantly X2X2' while Point D has about 50% X2X2 an d 35% X1X2 annihilation. 

Figure 1 shows the resulting pt distributions for leptons from decays to singlinos for all four 
benchmark points. Clearly, cuts on lepton transverse momentum of even 10 GeV will remove the wast 
majority of events for points B-D. However, one should notice that the distributions have considerable 
tails beyond the simple mass difference Am, due to the boost of the x\- 

2.2 Monte Carlo analysis 

We perform a Monte Carlo simulation of the benchmarks described above by generating both SUSY 
signal and SM background events with PYTHIA 6.413 [23]. The generated events are then put through 
a fast simulation of a generic LHC detector, AcerDET-1 . [24]. Although PYTHIA does not contain 
a framework for generating NMSSM events per se, it has the capability to handle the NMSSM spectrum 
and its decays. Since our scenario predicts the same dominant cross section as in the MSSM, namely 
gluino and squark pair-production, with negligible interference from the non-minimal sector, we use the 
built-in MSSM machinery for the hard process, and generate only squark and gluino pair-production. 

The detector simulation is done with standard AcerDET settings, with one exception: for detect- 
ing decays to the singlino the detector response to soft leptons is vital. We therefore parametrise the 
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efficiency of electron and muon identification as a function of lepton pr- For muons we base ourselves 
on the efficiencies shown in Figure 8-5 and 8-9 of the ATLAS TDR [25]; for electrons we use the same 
parametrisation scaled down by 0.82. While this is certainly not a perfect description of the real AT- 
LAS or CMS efficiencies during data taking, it incorporates some of the most important effects in an 
analysis, such as an absolute lower limit for lepton identification, at around 2-3 GeV for muons, and a 
difference in electron and muon efficiencies. However, it does not address other important issues, e.g. 
mis-identification of charged pions as electrons. To improve on these simple assumptions one would 
need a full simulation of the detectors, or efficiencies from data, which is clearly beyond the scope of 
this contribution. 

We generate events corresponding to 10 fb _1 of both signal and background (some with weights). 
Our background consists of large p T -binned samples of QCD 2^2 (10M), M^+jet (4M), Z+jet (3M), 
WW/WZ/ZZ (1M) and ti (5M) events. For the signal, PYTHIA gives a LO cross section of 24 pb for 
squark and gluino pair-production, thus 240 000 events are generated per benchmark point. 

We begin our analysis along the lines of the 'standard' di-lepton edge analysis. To isolate the 
SUSY signal from SM background we apply the following cuts: 

• Require at least three jets with pt > 150, 100, 50 GeV. 

• Require missing transverse energy $t > max(100 GeV, 0.2M e g ), where the effective mass M e g 
is the sum of the px of the three hardest jets plus missing energy. 

• Require two OSSF leptons with p T > 20, 10 GeV. 

After these cuts the background is small and consists mainly of ti, with some vector boson events sur- 
viving. The resulting di-lepton invariant mass distributions for points B and D can be seen in the left and 
right panels, respectively, of Fig. 2. The contribution from decays of to x\ via right and left-handed 
sleptons are shown in red and blue, other SUSY events, where the leptons mainly come from chargino 
or stau decays, in light grey, and the remaining SM background in dark grey. For the B benchmark point 
there is a small excess of events coming from decays of \% to singlinos (yellow) at low invariant masses, 
that survives due to other harder leptons in the event. However, all such events are removed for bench- 
mark point D because of the hard lepton px cut. 2 In this case one would miss the singlino and take the 
%2 to be the LSP dark matter candidate. 

It is clear that to increase sensitivity to the disguised NMSSM scenario, one needs to lower the 
lepton pt cuts. However, this opens the possibility for large increases in background. While most of this 
background, from uncorrected leptons, can in principle be removed by subtracting the corresponding 
opposite sign opposite-flavour (OSOF) distribution, assuming lepton universality, large backgrounds will 
increase the statistical error and a soft lepton sample is more vulnerable to non-universality from e.g. pion 
decays. The result of completely removing the pt requirement on the leptons is shown for benchmark 
points B and D in the left and right panels of Fig. 3, respectively. While there is indeed an increase in 
backgrounds, the effect on the signal is much more significant. For both benchmarks, the decay to the 
singlino is now visible as a large excess at low invariant masses. We have also tested scenarios with 
smaller values of Am, and find that we have a significant excess down to Am ~ 0.6 GeV, with the 
assumptions on lepton efficiencies described above. 3 

In the standard di-lepton analysis the edges of the red and blue distributions shown in Fig. 2 can be 
used to determine the relationship between the neutralino and slepton masses, in our scenario m? — fn~ 

and m~ — m~ . We extract additional information by also determining the position of the edge at low 
invariant masses, fitting a Gaussian-smeared step function to the OSOF subtracted distribution, shown 
in Fig. 4. In subtracting the OSOF distribution we have taken into account the asymmetry induced by 
the the difference in electron and muon efficiencies. This fit determines the mass difference Am, since 

2 Green denotes lepton combinations with one lepton from a slepton decay chain and the other from a decay to a singlino. 
The few pure singlino events at higher invariant masses are due to mis-combinations of leptons from different X2 decays. 
3 In fact, for such small mass differences we may also see displaced vertices due to the long lifetime of the x°- 
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m„ [GeV] m,, [GeV] 

Fig. 2: Di-lepton invariant mass distributions for point B (left) and point D (right) with standard lepton pr cuts. See text for 
colour coding. 




Fig. 3: Di-lepton invariant mass distributions for point B (left) and point D (right) without lepton pr cuts. 
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m™ ax = m^o — m^o. For point B the result of the fit is m™ ax = 2.93 ± 0.01 GeV, to be compared with 
the nominal value of 3.05 GeV, while for point D the result is m]^ ax = 0.77 ± 0.02 GeV, with a nominal 
value of 0.87 GeV Both results are significantly on the low side with respect to the small statistical 
errors. We speculate that this systematic error is at least in part due to the step function used in the fit to 
the edge, and that a more sophisticated description will give results closer to the nominal values. 

A final comment is in order concerning early discoveries: in fact, since all SUSY cascades will 
contain the decay to a singlino, the lower edge in the di-lepton distribution may appear much earlier than 
the 'standard' decay through a slepton, if at all present, provided that the soft leptons are searched for. 



2.3 Conclusions 

We have demonstrated that lowering the requirements on lepton transverse momentum in the standard di- 
lepton edge of SUSY searches may reveal unexpected features, such as the NMSSM in disguise. While 
our numerical results are sensitive to the exact lepton efficiencies, to be measured at the experiments, 
and while there may be additional backgrounds not simulated, such as multi gauge boson and/or multi jet 
final states, the OSOF subtraction procedure ensures that the background is removable and the NMSSM 
scenario in question is discoverable down to very small mass differences Am = m^o — m^o . 
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3. The MSSM with decoupled scalars at the LHC 1 

3.1 Introduction 

Assuming a large soft-breaking scale for all MSSM scalars [26, 27, 28, 29, 30] pushes squarks, sfermions 
and heavy Higgses out of the reach of the LHC without affecting the gaugino sector. Even though the hi- 
erarchy problem will not be solved without an additional logarithmic fine tuning in the Higgs sector, such 
models can be constructed to provide a good dark-matter candidate and realize grand unification while 
minimizing proton decay and FCNCs. We investigate their LHC phenomenology, with all scalars decou- 
pled from the low-energy spectrum. We focus on gaugino-related signatures to estimate the accuracy 
with which its underlying parameters can be determined [31, 32]. 

3.2 Phenomenology 

The spectrum at LHC mass scales is reduced to the Standard-Model with a light Higgs, plus gauginos 
and Higgsinos. At the high scale Mg the effective theory is matched to the full MSSM and the usual 
renormalization group equations apply. The Higgsino mass parameter fi and the ratio tan (3 in the Higgs 
sector correspond to their MSSM counter parts. The gaugino masses Mi 2,3 and the Higgs-sfermion- 
sfermion couplings unify, and Mg replaces the sfermion and the heavy Higgs' mass parameters. This set 
resembles the mSUGRA parameter set except for tan [3 now playing the role of a matching parameter 
(with all heavy Higgses being decoupled) rather than that of an actual vev ratio [33]. 

We select our parameter point lead by three constraints: first, we minimize the amount of fine 
tuning necessary to bring the light Higgs mass into the 100 to 200 GeV range and reduce Ms to 10 TeV, 
which is still outside the LHC mass range. Another reason for this low breaking scale is that we want the 
gluino to decay inside the detector (preferably at the interaction point) instead of being long-lived [34, 
29]. Heavier sfermions increase the life time of the gluino such that it creates a displaced vertex or even 
hadronizes [35]. 

Secondly, we obtain the correct relic dark-matter density VLh 2 = O.lll^Q QQg [36] by setting 
H = 290 GeV and M 2 (M GU t) = 132.4 GeV or M 2 (M wcak ) = 129 GeV. This corresponds to the 
light-Higgs funnel ttilsp ~ M 2 /2 « m^/2, where the s-channel Higgs exchange enhances the LSP 
annihilation rate. And finally, rrih needs to be well above the LEP limit, which we achieve by choosing 
tan P = 30. We arrive at a parameter point with = 129 GeV, m g = 438 GeV, chargino masses 
of 117 and 313 GeV, and neutralino masses of 60, 117, 296, and 310 GeV, using a modified version of 
SuSpect which decouples the heavy scalars from the MSSM RGEs [30, 37]. The neutralinos/charginos 
X% an d X\ as wei l as Xa an d X2 aTe degenerate in mass. All neutralinos/charginos and most notably the 
gluino are much lighter than in the SPSla parameter point, which greatly increases all LHC production 
cross sections. It is important to note that this feature is specific to our choice of parameters and not 
generic in heavy-scalar models. 

Table 1 shows the main (NLO) cross sections at the LHC [38, 39]. The SUSY production is 
dominated by gluino pairs whose rate is eight times that of the SPSla point: the lower gluino mass 
enlarges the available phase space, while in addition the destructive interference between s and i-channel 
diagrams is absent. The second largest process is the xfx® production, which gives rise to a 145 fb of 

'R. Lafaye, T. Plehn, M. Rauch, E. Turlay, and D. Zerwas 



99 


63 pb 


x^g 


0.311 pb 


X ± X° 


12 pb 


X°~9 


0.223 pb 


X ± X T 


6pb 


x°x° 


0.098 pb 


Total 






82 pb 



Table 1 : NLO cross sections for SUSY pair production at the LHC. Branching ratios are not included. 
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observables 


systematic error 


statistical error 


theory error 


m h 


128.8 GeV 


0.1% 


energy scale 


0.1% 


4% 


X.2 Al 


57 GeV 


0.1% 


energy scale 


0.3% 


1% 




145.2 fb 


5% 


luminosity 


3% 


20% 


g — > 6/not(6) 


0.11 


5% 


b tagging 


0.3% 


20% 




68.2 pb 


5% 


luminosity 


0.1% 


20% 



Table 2: Summary of all observables and their errors. We assume an integrated luminosity of 100 fb 



hard-jet free, e and [i trilepton signal, more than a hundred times that of the SPSla parameter choice. 
3.3 Observables 

The first obvious observable is the light Higgs mass m^. Although slightly higher than in most MSSM 
points, rrih can still be measured in the Higgs decay to two photons [40]. The systematic error on this 
measurement is mainly due to the uncertainty of the electromagnetic energy scale. 

A measurement of the gluino pair production cross section appears feasible and can be very helpful 
to determine M3 [41]. Most gluinos (85%) will decay through a virtual squark into a chargino or a 
neutralino along with two jets. The chargino will in turn decay mostly into the LSP plus two leptons or 
jets. Such events would then feature at least 4 hard jets, a large amount of missing energy and possibly 
leptons. The main backgrounds are ti pairs (590 pb) and VF+jets (4640 pb) as well as Z+jets (220 pb). 
Despite these large cross sections, we have checked using a fast LHC-like simulation that most of the 
background can be eliminated by requiring a minimal number of hard jets or by applying standard cuts 
on the missing energy or the effective mass M e g =I$t + J^PTj- The main source of systematic errors 
on the cross section is the 5% error on the luminosity. The theory error on the cross section we estimate 
to 20%. 

The next relevant observable is the trilepton signal. After gluino pairs, the second-largest rate 
comes from the direct production of xfx% Wltn 22% of xf s decaying through a virtual W into an elec- 
tron or muon, a neutrino and the LSP. Similarly, 7% of decay through a virtual Z into an opposite- 
sign-same-flavor lepton pair (OSSF) and the LSP The resulting signal features three leptons (two of 
them with OSSF), missing energy from the LSPs and the neutrino, and no jet from the hard process. The 
backgrounds are mainly WZ (386 fb) and ZZ (73 fb), the latter with one missed lepton. Taking into 
account all branching ratios [42], the trilepton signal has a rate of 145 fb. Without any cuts, the identifi- 
cation efficiencies of 65% (e) and 80% (fi) leave us with 1 10 to 21 1 fb for the background and 40 to 74 fb 
for the signal, depending on the number of electrons and muons in the final state. A dedicated study with 
the appropriate tools would evidently provide a better understanding of signal and background. As in the 
previous case, the main source of systematic errors is the luminosity. We also take the theory error on 
the value of the trilepton cross section to be roughly 20%. 

For the trilepton signal we can define a kinematic observable: 10% of X2 S decay into an OSSF 
lepton pair and the LSP. The distribution of the invariant mass of the leptons features a kinematic upper 
edge at m^o — m^o. Such an observable gives precious information on the neutralino sector and hence 
on Mi. Its systematic error is dominated by the lepton energy scale. The statistical error we estimate to 
be of the order of 1%, from a ROOT fit of the Mu distribution. Finally, we use the ratio of gluino decays 
including a b quark to those not including a b. We roughly assume a systematic error of 5% due to the b 
tagging and 20% on the theory prediction. 

Table 2 summarizes the central values and errors for all observables we use. The third and fourth 
columns give the experimental systematic errors and their sources, the fifth column gives the statistical 
errors corresponding to a few years of the LHC's nominal luminosity (100 fb -1 ). The last column gives 
a conservative estimate of the theory uncertainties. 
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3.4 Parameter determination 

To study the effects of the different error sources, we first look at a low-statistics scenario and ignore 
all theory uncertainties. Then, we choose the limit of high statistics to estimate the ultimate precision 
barrier imposed by experimental systematical errors. Finally, we look at the effect of theory errors by 
including them into the previous set. We expect the theory errors to dominate, based on the currently 
available higher-order calculations. We use the parameter extraction tool SFitter [31], which in parallel 
to Fittino [32] was developed for this kind of problem. 

With no information on the squark and sfermion sector, except for non-observation, we are forced 
to fix Ms and A t in our fits. Moreover, we set Mj = Mi at the unification scale Mgijt, lacking enough 
information from the neutralino/chargino sector. Using Minuit, we then fit the remaining parameters to 
the LHC observables. The x 2 minimum we identify using Migrad, while Minos determines the approxi- 
mately Gaussian errors. Note that in a more careful SFitter study we would use flat theory errors [43, 31], 
but given the huge difference in computing time we employ a Gaussian approximation in this preliminary 
study. Our distant starting point is (Mi, M 3 , tan (3, fi) = (100, 200, 10, 320). 

Table 3 shows the result of the different fits. It is interesting to note that tan (3 is undetermined 
except in the case of infinite statistical and theory's accuracy. This is due to the fact that only one of 
the five Higgs masses is measured. We present a study on the determination of tan (5 from (g — 2) M 
elsewhere in this volume. The quality of the trilepton and gluino signals gives very good precision on 
the determination of Mi and M3, even with low statistics. Including theory errors indeed decreases the 
accuracy but still allows for a determination of the mass parameters: M3 only depends on the large gluino 
rate and its decays, explaining its relative stability for smaller statistics. The weakly interacting Mi and 
M2 suffer a larger impact from the theory errors, because they depend on the trilepton rate and also on 
the b to non-6 gluino-decay ratio, both of which bear a large theory error. 

3.5 Outlook 

The MSSM with heavy scalars is built to satisfy current experimental and theoretical constraints on 
physics beyond the Standard Model while keeping some of the features of the TeV-scale MSSM. At the 
LHC, light gauginos and Higgsinos will lead to sizeable production rates, allowing us to study these new 
states. 

The main observable channels are gluino pairs and the tri-leptons, whose hard-jet free channel 
makes it a fairly clean channel with respect to Standard-Model and SUSY backgrounds. Additional 
observables such as the light Higgs mass, the {m^o — m^o) kinematic edge and the 6-to-non-fe gluino- 
decay ratio give us access to most parameters at the level of a few percent with 100 fb _1 luminosity, 
based on experimental uncertainties. Theory errors increase the error bands on the model parameters to 
0(15%). 

Obviously, the scalar sector including tan (3 is only poorly constrained, if at all. New comple- 
mentary observables could improve this limitation. Similarly, a look at other parameter points would be 



parameter 


nominal 


fitted 


low stat. 


00 stat. 


00 stat.+theory 


M 2 


132.4 GeV 


132.8 GeV 


6 


5% 


0.24 


0.2% 


21.2 


16% 


M 3 


132.4 GeV 


132.7 GeV 


0.8 


0.6% 


0.16 


0.1% 


5.1 


4% 




290 GeV 


288 GeV 


3.8 


1.3% 


1.1 


0.4% 


48 
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tan/3 
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28.3 
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1.24 
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132.4 GeV 


132.8 GeV 
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fixed 






M S 


lOTeV 






fixed 







Table 3: Result of the fits. Errors on the determination of the parameter are given for the three error sets described in the text. 
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needed to remove the specific properties of the point studied and to provide a more complete view of the 
LHC discovery potential of a MSSM with decoupled scalars. 

Acknowledgements 

We would like to thank the organizers of the great Les Houches 2007 workshop Physics at TeV-Scale 
Colliders. 



23 



4. Finding the SUSY mass scale 1 
4.1 Introduction 

In the past it has been suggested that a good starting point for the determination of the mass scale of new 
susy or other exotic particles is the "effective mass" distribution [44, 45]. There are a number of slightly 
different definitions of Mes and the phrase "mass scale" [46] but a typical definition of Mes would be 

M ES = p ™° + J2pT(i), (1) 

i 

in which p™ lss is the magnitude of the event's missing transverse momentum and where Px(i) is the 
magnitude of the transverse momentum of the i-th hardest jet or lepton in the event. 

All definitions of Mes are motivated by the fact that new TeV-scale massive particles are likely 
to be produced near threshold, and so by attempting to sum up the visible energy in each event, one can 
hope to obtain an estimate of the energy required to form the two such particles. Broadly speaking, the 
peak in the Mes distribution is regarded as the mass-scale estimator. 

Although the effective mass is a useful variable, and simple to compute, it has undesirable proper- 
ties: The desired correlation between Mes and the mass scale relies on the assumption that the particles 
are produced near threshold. While it is true that the cross sections will usually peak at threshold, they 
can have significant tails extending to y/i values considerably beyond the threshold value. It is very hard 
to make precise statements about the mass scale from Mes alone. 

In this letter we introduce a variable, rriTGen, which is designed to make more precise measure- 
ments of the mass scale by using event kinematics, rather than simple energy sums or ad-hoc rules. The 
aim is to produce event-by-event lower bound on the mass of pair-produced heavy particles. The variable 
has been constructed so that it is our ignorance about the ancestry of the final state particles, and the loss 
of information from the invisible massive heavy particles. 

A solution already exists for the simplest case of interest - in which the final state consists of only 
two visible particles (plus the invisibles). For this case there is no combinatorial problem, as one can 
assume that each visible particle belongs to one of the initial heavy particles (e.g. one jet comes from 
each squark parent). The variable defined for that case is known as rriT2 and is described in [47, 48]. 

The generalisation to the case of arbitrary numbers of final state visible particles is the subject 
of this letter and is called rriTGen- rriTGen is defined to be the smallest value of obtained over 
all possible partitions of momenta in F into two subsets a and (3 - each subset representing the decay 
products of a particular "side" of the event. Note that mn is itself defined in terms of p T and m a 
(respectively the transverse momentum and mass of one side of the event), and mp (respectively the 
transverse momentum and mass of the other side of the event), and \ (the mass of each of the unobserved 
particles which are supposed to have been produced on each side of the event) as follows: 

m T 2 (Pt , Pt > Pr 1SS ,m a ,rrip,x) = 

maxjm^p^q^jm^x), m T (p T , q^ 2) ; m p , x)}] • (2) 



mm 



where 
in which 



m 2 T (PT, p£ ; m«, X) = m 2 a + X 2 + 2{E%E* - p T ■ P T ) (3) 



£t = V(Pt) 2 + < and ^ = V(Pt) 2 + X 2 (4) 

and likewise for a < — ► j3. With the above definition (in the case x = m x°)> m T2 generates and event- 
by-event lower bound on the mass of the particle whose decay products made up either of the two sides 
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of the event, under the assumption that the event was an event represents pair production followed by 
decay to the visible particles and an unseen massive particle on each side. When evaluated at other values 
of x the above properties are retained approximately (see [47, 48, 49, 50, 51, 52]). There exist events 
which allow this lower bound to saturate, and so (in the absence of background) the upper endpoint of 
the ?7iT2 distribution may be used to determine the mass of the particle being pair produced. 

4.2 Example distributions 

In this letter we show the results for simulations of an example particle spectrum for proton-proton col- 
lisions at LHC centre-of-mass energy of yfs =14 TeV. The HERWIG [53, 54, 55] Monte Carlo generator 
was used to produce inclusive unweighted supersymmetric particle pair production events. Final state 
particles (other than the invisible neutrinos and neutralinos) were then clustered into jets by the longi- 
tudinally invariant clustering algorithm for hadron-hadron collisions [56] used in the inclusive mode 
with R = 1.0 [57]. Those resultant jets which had both pseudo-rapidity (rj = — In tan 9/2) satisfying 
\t]\ < 2 and transverse momentum greater than 10 GeV/c were used to calculate rriTGen and M^g. 




Fig. 1: On the left hand side is a graphical representation of the susy mass spectrum of the point described in the text. The 
vertical positions of the particles indicate their masses. The horizontal positions of the centres of the bars indicate the relative 
LHC production cross-section (arbitrary units). The lines joining particles indicate decays with branching fractions in the 
following ranges: greater than 10 -1 solid; 10 -2 — > 10 -1 dashed; 10 -3 — + 10 -2 dotted. The middle plot shows the distribution 
of our variable, MrGen, with MrGen increasing vertically to ease comparison with the spectrum. The right hand plot shows 
the distribution of another variable, Mes/2, where Mej is defined in eq. 1. In both the MrGen and the Mej plots, the lighter 
shading shows the histograms with the number of events multiplied by a factor of twenty, so that the detail in the upper tail may 
be seen. 

In figure 1 we show the distributions of Mxcen and M^s for a sample point with a spectrum, 
defined by the mSUGRA parameters: {m = 1200 GeV, mi = 420 GeV, tan/3 = 10, m t = 174 GeV, 

2 

(i < 0}, The spectrum and branching ratios were calculated using Isa jet [58] version 7.58. In these 
plots we assume it is possible to accurately assign all visible momenta to the correct category F or G, 
i.e. "interesting final state momenta" versus "initial state radiation". Plots in which this is not assumed 
look similar and can be found in [59]. The HERWIG initial state radiation and underlying event have 
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(a) (b) 



Fig. 2: A/tGeu distributions for a Les Houches blind data sample as described in the text, (a) Invisible particle mass, x — 0. 
(b) The colour scale shows the number of events with a particular value of Mroen (y-axis) as a function of the input invisible 
particle mass, \ (x-axis). 

been switched off, and the parameter x which is required to calculate Mxcen has been set to the mass of 
the lightest supersymmetric particle. 

It can be seen that, as intended, the upper edge of the distributions gives a very good indication of 
the mass of the heaviest pair-produced sparticle. Other supersymmetric points show similar behaviour 
[59]. This means that the position of the upper edge of rriTGen can be used to find out about the mass 
scale of any semi-invisibly decaying, heavy, pair-produced particles. 

Furthermore, a change in slope can be observed at lower masses due to significant pair production 
of lower-mass particles (e.g. chargino and/or neutralino pairs). Therefore it is also possible in principle 
to extract from information about several different mass scales. 

In figure 2 we plot the rriTGen distribution for events from the Les Houches 2007 blind data 
sample 2 . We select events which have at least one jet with |^| < 3.2 and transverse momentum greater 
than 400 GeV, and also have missing transverse momentum greater than 70 GeV. We veto any event 
containing an electron, muon or photon with transverse momentum greater than 20 GeV. 3 There is no 
evidence for an end-point in the distribution, like the one seen in figure 1 or in the other examples in 
[59]. Neither is there evidence for any "kinks" in the rriTGen distribution when plotted as a function of 
the invisible particle mass, \ (figure 2b). Such end-points and kinks would be expected if the events 
contained many-body or cascade decays of strongly interacting objects to visible and invisible heavy 
particles [49, 50, 51, 52], and so disfavours a model of this type. 

4.3 Conclusions 

In conclusion, rriTGen could be a useful variable for determining mass scales at the LHC. If you are 
interested in using it, contact the authors for the code and let us know what you find. 
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5. A hybrid method for SUSY masses from fully identified cascade decays 1 

5.1 Introduction 

This letter describes a new hybrid technique for improving the precision with which SUSY particle 
masses can be measured at the LHC. Existing techniques usually make use of the positions of end- 
points in experiment-wise distributions of invariant mass combinations of visible SUSY decay products 
[60, 61, 62], or use E™ tss constraints from 'symmetric' events in which the same SUSY decay chain 
has participated in both 'legs' of the event [48]. In both cases some information regarding the events 
is discarded - in the former case event-wise E™ lss information is not used in the experiment-wise end- 
point analysis, in the latter experiment-wise invariant mass end-point constraints are not used in the 
event-wise analysis. In this letter we describe a simple 'hybrid' technique which enables optimum use of 
both experiment-wise and event-wise information to fully reconstruct SUSY events and hence improve 
the mass measurement precision. 

5.2 Description of technique 

The new technique involves conducting a kinematic fit to each selected SUSY event, with the sparticle 
masses as free parameters. Crucially, the x 2 function of the fit involves both event-wise E™ lss constraints 
and experiment-wise invariant mass end-point constraints. It should be appreciated that without the 
Ejniss constraints each event-wise fit is formally equivalent to solving the experiment-wise invariant 
mass end-point constraints for the individual sparticle masses, and consequently each fit will give the 
same value for each mass. The RMS values of the distributions of these masses will then be consistent 
with the mass precisions obtained from the conventional method involving solution of the end-point 
constraints. Addition of the event-wise E™ lss constraints reduces the number of degrees-of-freedom of 
the kinematic fits and hence can improve the mass measurement precision. In this case the widths of the 
distributions of mass values obtained from different events for one Monte Carlo experiment are larger 
than those obtained when E™ tss constraints are excluded, however the means of the distributions across 
many such experiments measure the masses more accurately. 

5.3 Example: qi decays in SPSla 

At mSUGRA point SPSla there is a significant branching ratio for the decay chain 

QL - X° 2 q - l R lq - xlllq. (D 

This chain provides 5 kinematic end-point mass constraints from invariant mass combinations of jets and 
leptons [63]: 

• m(ll) max = 77.08 ± 0.08(scale) ± 0.05(stat) GeV 

• m(llq) max = 431.1 ± 4.3(scale) ± 2.4(stat) GeV 

• m{llq) min = 203.0 ± 2.0(scale) ± 2.8(stat) GeV 

• m(lq)™ x = 380.3 ± 3.8(scale) ± 1.8(stat) GeV 

• m(lq)™ ax = 302.1 ± 3.0(scale) ± 1.5(stat) GeV 

For this study unbiased samples equivalent to 100 fb _1 (one Monte Carlo 'experiment') of SPSla 
signal events and ti background events were generated with HERWIG 6.4 [64, 54] and passed to a 
generic LHC detector simulation [65]. A lepton reconstruction efficiency of 90% was assumed. 

Events were selected in which the above decay chain appears in both legs of the event with the 
following requirements: 

• N jet > 2, withp T (j2) > 100 GeV, 
'MM. Nojiri, G. Polesello and D.R. Tovey 
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• M e//2 = E™ ss + PT (jl) +Pt(j'2) > 100 GeV, 

• EJp iss > max(100GeV,0.2M e//2 ), 

• Ni ep = 4, where lep = e//x(isolated) and pt(W) > 6 GeV, 

• 2 Opposite Sign Same Flavour (OSSF) lepton pairs. If the pairs are of different flavour both pairs 
must have m(U) < m(ll) max . If both pairs are of the same flavour then one and only one of 
the two possible pairings must give two m(ll) values which are both less than m(ll) max . These 
pairings allocate the leptons to each leg of the event. 

• One and only one possible pairing of the two leading jets with the two OSSF lepton pairs must 
give two m(llq) values less than m{llq) max . These pairings allocate the jets to each leg of the 
event. 

• For each inferred leg of the event the maximum(minimum) of the two m(lq) values must be less 
than m(lq)™^ x o y This ordering allocates the leptons to the near and far [61] positions in the 
decay chain. 

The requirement of 4-leptons in two OSSF pairs and two high-p-p jets consistent with kinematic 
end-points, together with large E™ lss , is effective at removing the majority of SM and SUSY back- 
grounds (see below). 

Each selected event was fitted with MINUIT [66]. Free parameters were taken to be the four 
masses appearing in the decay chain: m(qx), ^(x^)' m (^) an ^ The mass-shell conditions 

and measured momenta of the visible decay products for each leg were solved to determine the LSP 
four-momenta, giving two solutions for each leg. The \ 2 minimisation function was defined by: 

^ 2 = / m(ll)? v r-m(li)™ x t " 



Vm{liy 

2 



+ . m (%g-^g y + /«-^sV 



+ 



ax x 
;expt 



'rn(lq)r? vt -m(lq)r*\ 2 ^ ( m(lq)™ x vt - m(lq)r x 



a m{lq)™°- x I \ a m(h)To a 





p x {x\{i))+Px{x^))-Er ss \ , (py(m))+Py(m))-E- 



+ fMAn " ' + — " — " »— , (2) 



(7 JTJmiss J y (7 J^miss 



where evt denotes an expected end-point value derived from the masses in the event-wise fit with the 
formulae of Ref. [61], and expt denotes a 'measured' experiment- wise end-point value. The uncertain- 
ties a in these 'measured' endpoints were those quoted above. The uncertainties on the measurements 
of the x and y components of EJp iss , o E miss and a E rniss, were given by 0.5y / £ i |! 7 " r where E^ m is the 
scalar sum of jet p? of the event. This function incorporating both event -wise E?p l ss constraints and 
experiment-wise end-point constraints was evaluated for each of the four pairs of Xi momentum solu- 
tions obtained from solving the leg mass-shell conditions. Fitted masses were obtained when x 2 was 
minimised for the event. Fitted masses were used in the subsequent analysis only if MINUIT judged the 
fit to have converged and Xmin < 35.0. 

Following application of the selection cuts described above and the requirements of fit convergence 
and low fit Xmin 38 SUSY 'signal' events with the above decay chain appearing in both legs were 
observed. 4 SUSY background events were observed, consisting of the above decay chain in both legs 
but with one or two leptonically decaying staus produced in the decays of the x^'s- No ti background 
events were observed in 100 fb _1 equivalent data. More SM background events may be expected in a 
real experiment, given that effects such as charge and lepton mis-identification are not included in the 
fast detector simulation. Nevertheless the contribution is still expected to be negligible. 
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Fig. 1 : Distributions of sparticle masses. The four distributions on the left are of masses obtained from event-wise fits, for one 
MC experiment. Each entry is obtained by minimising the \ 2 function shown in Eqn. 2. The four distributions on the right are 
likelihood distributions of sparticle masses obtained from 100 MC experiments. Each entry is the mean of an experiment-wise 
mass histogram such as those on the left. 



Each event-wise fit generated one set of values for the sparticle masses, namely those values which 
minimise the broad \ 2 function in Eqn. 2. The distributions of these values for one Monte Carlo exper- 
iment are shown in Fig. l(left). In order to demonstrate the performance of the technique and judge the 
uncertainties in the measurements the above procedure was repeated for 100 Monte Carlo experiments. 
For each experiment, kinematic end-point positions were sampled from gaussians with means and sig- 
mas given by the means and uncertainties listed above. The five sampled end-point positions for each 
experiment were solved simultaneously with a MINUIT fit to give initial mass values for input to the 
MINUIT event-wise kinematic fits. For each experiment relative jet(lepton) energy scale values were 
sampled from gaussians of width 1%(0.1%) reflecting likely ultimate energy scale uncertainties at the 
LHC. Each experiment generated a set of sparticle mass histograms similar to those shown in Fig. l(left). 
The means of these histograms for the 100 MC experiments were then used to construct likelihood his- 
tograms for the masses, shown in Fig. l(right). The standard deviations of these histograms were taken 
to provide the uncertainties on the sparticle mass measurements. 

Unbiased MC data equivalent to only one 100 fb _1 experiment were available for this study. 
For this reason the same events were used for each MC experiment, with just the end-point values and 
jet/lepton energy scales varying. The additional uncertainties in the final mass values expected from 
varying event samples were estimated from the mean statistical uncertainties in the mean experiment 
mass values as extracted from the event-wise distributions such as those shown in Fig. l(left). We 
evaluated the experiment-by-experiment spread due to varying event samples as a/y/n, where a is the 
RMS of the event-wise distributions as shown in Fig. l(left), and n is the number of entries in each plot. 
These additional contributions were added in quadrature to the uncertainties obtained from the study. 
This approximation was checked with a second sample of SPSla events equivalent to 100 different MC 
experiments, biased to force gluinos to decay to qi, b or t, qi to decay to \\ an( ^ X% to decay to e or fl. 

The results of this study are summarised in Table 1. For comparison purposes the analysis was ini- 
tially carried out with the E™ %ss constraints removed from the \ 2 function. The measurement precisions 
are consistent with those obtained from the conventional end-point fitting method, as expected following 
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State 


Input 


End-Point Fit 


Hybrid Method, E™ ss 


Hybrid Method, no E™ lss 






Mean 


Error 


Mean 


Error 


Mean 


Error 




96.05 


96.5 


8.0 


95.8(92.2) 


5.3(5.5) 


97.7(96.9) 


7.6(8.0) 




142.97 


143.3 


7.9 


142.2(138.7) 


5.4(5.6) 


144.5(143.8) 


7.8(8.1) 


v° 
A2 


176.81 


177.2 


7.7 


176.4(172.8) 


5.3(5.4) 


178.4(177.6) 


7.6(7.9) 


<7L 


537.2-543.0 


540.4 


12.6 


540.7(534.8) 


8.5(8.7) 


542.9(541.4) 


12.2(12.7) 



Table 1: Summary of mass measurement precisions for SPSla states. Column 2 lists masses used in the HERWIG generator, 
Columns 3 and 4 the fitted masses and uncertainties obtained from the conventional fit to kinematic end-points, Columns 5 
and 6 the equivalent values obtained with the new technique and Columns 7 and 8 the equivalent values obtained with the 
new technique excluding E™ 133 constraints. Figures in parentheses are those obtained with the biased sample of non-repeated 
events. All masses are in GeV. The quoted mass range for qL excludes b squarks, which are produced less readily than the light 
squarks. 



the reasoning outlined above. The analysis was then repeated including the E™ zss constraints, giving 
an overall improvement in sparticle mass precisions ~ 30% for all four masses considered. A similar 
improvement was found when using the biased sample of non-repeated events for different experiments. 
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6. A blind SUSY search at the LHC 1 

6.1 Introduction 

Most of the SUSY prospects [67] [68] at the LHC are based on simulations and analyses where the 
analyzer knows all the properties of the signal as well as those of the SUSY and Higgs backgrounds to 
this signal. 

In constrast with these situations, we describe in this letter a new SUSY data challenge at the LHC which 
is called the "Blind SUSY Search Project" [69]. As name of this project suggests, here, the analyzer 
ignores the properties of the searched SUSY signal. 

We have produced a 100 pb~ l pseudo-data (PsD{) sample which consists in a randomized mixture of 
the Standard Model (SM) backgrounds and the inclusive Higgs and SUSY production of an unrevealed 
SUSY model. We also provide separate and independent samples of the SM backgrounds. 

The aim of the challenge is to determine the type of underlying SUSY breaking mechanism as well as 
the corresponding parameters of the hidden SUSY model. 

All of the samples are under a simple ROOT [70] format so as to propose this challenge not only to the 
experimental HEP community but also to the theorists. 

The motivations for this challenge lie in the following questions. Let's hypothesize the presence of a 
SUSY signal at the LHC, can one: 

• determine the excess with respect to the SM expectations and quantify it? 

• handle the possible presence of several SUSY and Higgs signals and how does that affect the measure- 
ment of experimental observables (masses, mass differences, cross sections, ratios of branching fractions, 
spins,...)? 

• determine the type of underlying SUSY breaking mechanism at play (gravity or gauge or anomaly 
mediated) 

• distinguish different types of phenomenological hypotheses (R-parity conservation, phases, high scale 
unifications,...)? 

• evaluate the values of the parameters of the underlying SUSY model? 

Part of such questions has been posed and partially answered in previous SUSY challenges. But most 
of them provided either an exclusive signal (ie: without all the decay channels open), or were missing 
the Higgs and SUSY background, or part of the SM backgrounds,... In the current challenge we tried to 
provide all of these pieces. 

In section 2, we'll describe in some details how the samples were produced. In sections 3 and 4 respec- 
tively we'll explain how to access and how to analyze the data. 

6.2 Samples production 

Production tools 

All the processes were generated using Pythia version 6.325 (v6.325) [71]. The SUSY mass spectrum 
and decay table in the vl.O "SUSY Les Houches Accord" (SLHA)[72] format was read in by Pythia. 
The CTEQ6L1 [73] proton parton density functions (PDF) were utilized for all the processes through an 
interface to the LHAPDF [74] v5.2.3 package. 

For the SM backgrounds and the Higgs processes the r were decayed by Tauola [75] v2.6. However, 
because of missing pieces in the Tauola interface, the r from SUSY processes were decayed by Pythia 
unabling to account for the spin correlations in these cases. 

The ATLAS detector response was simulated using a personal fast simulation based on ATLFAST [76] 
v00-02-22. All the reconstructed quantities are simulated using smearing functions and the correspond- 
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ing variables are also used to mimic the trigger conditions 2 as described in the ATLAS High Level Trigger 
TDR [77]. Note that the events failing the trigger conditions were removed from the PsD\ set, but not 
from the background samples. Obviously the Monte Carlo (MC) truth informations were removed from 
the PsD\ set as well. 

The output of the fast simulation is a PAW [78] ntuple which is subsequently converted into a root-tuple 
using ROOT v5. 14. 

For the samples normalization an integrated luminosity of 100 pb^ 1 was assumed and the Pythia leading 
order (LO) cross sections were used for the sake of simplicity. Nevertheless some background (sub)- 
processes with very high cross sections had to be arbitrarily taken out of the PsD\ sample in order to 
keep the number of events to be produced within reasonable limits. Table 1 contains the full list of these 
pruned sub-processes: 



Process 


Pythia 
Process Index 


Kinematical Cut 


qq -> qq 


MSEL = 1 


Pt < 160 GeV 


(q = u/d/s/g) 






cc 


MSEL = 4 


p T < 40 GeV 


bb 


MSEL = 5 


p T < 40 GeV 


7 + jets 


M5C/ J B(14,29,115) = 1 


PT < 20 GeV 


Low mass resonances 






(*, T, *,...) 






Elastic Scattering 






Diffraction 







Table 1 : The list of high cross section processes that were removed from the PsDi sample. 



The PsD\ sample is stored in 7 root-tuples and contains in total 4.5M events for a total size of 12 Gb. 

Each background root-tuples contains exactly 100k events. In total there are 1593 such root-tuples for a 
total statistics of 159. 3M events that amount to 424 Gb. 

All the other details about the events generation can be found on the project website: [69]. 
6.3 Access to data 

The full dataset was too large to be kept on disk. Therefore it is stored on tape at the Lyon Computing 
Center, except for the PsD\ sample which can also be downloaded from the project website. 

For those who have an account at this Computing Facility, the samples are available on HPSS[79] in 

the cchpssdO /hpss/in2p3.fr/home/m/muanza/GDR_SUSY/SUSY_Blind/ directory. And the 
following sub-directories: Pseudo-DATA/ final /0 / 100 Jnvjpb/, BKGD/ < bkgdjprocess > and 
ANALYSIS/ SKIMMING /mETlhO / respectively contain the PsD^ the SM background and the 
skimmed (Et > 150 GeV) samples. 

For those who don't have an account at the Lyon Computing Center, there's a possible data access using 
SRB[80]. The participant has to be registered as an SRB user. This can be done by sending me an email 
at muanza@in2p3.fr. Then the participant needs to install an SRB client on his (her) computer to be 
able to list and copy the root-tuples. The main SRB directory is 

/ home /smuanza.ccin2p3/GDR_SUSY/ SUSY _Blind/ and all the sub-directories structure is that 
of HPSS. 

All the useful details about SRB are explained on the projcet website. 

2 The participants are asked to present results only with events passing at least one of these conditions. 
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6.4 Data analysis 

Template analysis program 

A template analysis tarball is provided on the project website 

(http : I j www — cluedO.f 'nal.gov/%7Emuanza/ ' BlindSU SY / Analysis/ Run -Analysis. tar. gz). 
It can either be used as is or be hacked by the participants. In any case, the most useful information for 
those who'd like to write their own analysis code are the cross sections of all the SM backgrounds. They 
can be found in the "proc-xsect" array at the top of the Analysis.C file. 

Some illustrations of the PsD\ sample 

Here are some plots advertising the project. They were produced after rejecting the events not passing 
any trigger requirements 3 . 




Fig. 1: The distributions of the missing (left) and scalar (right) Et after rejecting the events failing the trigger requirements. A 
SUSY signal (dark blue) studied in [81] and its SUSY background (light blue) are superimposed only for illustrative purposes. 

These plots exhibit an excess of pseudo-data events in the high missing and scalar Et tails... 
Signal templates production 

A tarball of the production package is available on the project website 

{http : //www — cluedO.fnal.gov/%7Emuanza/ BlindjSUSY/ Production/running^susyblind.tar.gz). 
The analyzers have to produce their SLHA input cards for their preferred SUSY models and to produce signal 
templates using the production tarball. This way they can test ideas about possible signals that could explain the 
difference between the PsD\ sample and the SM background. 

6.5 Conclusions and prospects 

We have proposed a special SUSY data challenge at the LHC that includes the full Higgs and SUSY inclusive 
production for an unrevealed SUSY model on top of the SM backgrounds. 

The aim of this challenge is to subject to a blind analysis possible sttategies to disentangle a given Higgs or SUSY 
signal in the presence of simultaneous contributions from different other Higgs and SUSY processes. And to see 
how well these strategies enable to determine the properties of the SUSY model under study. 

We look forward for participants to this challenge to present their analysis of this first pseudo-data sample. 
We are eager to see what experimental observables they'll have measured, what will be their uncertainties esti- 
mates. We are expecting their initial best guess for the values of the hidden SUSY model parameters. And we 
also suggest they provide the SUSY fitter groups [82] [83] with their observables and uncertainties so as to find out 
what global fits could teach us about in this blind analysis context. 



3 A code for simulating the trigger conditions can be found in the ana:HLT() function of the ana.C file in the analysis tarball 
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7. Off-shell effects for decay processes in the MSSM 1 

7.1 Introduction 

Theoretical arguments and experimental observations indicate that new particles or interactions play an important 
role at the TeV scale, which will become directly accessible at the Large Hadron Collider (LHC) and its planned 
complement, the International Linear Collider. In the near future we can therefore anticipate ground-breaking 
discoveries that reveal physics beyond the Standard Model (BSM) and allow to gain insight into the structure 
of the fundamental theory. Theoretically appealing extensions of the Standard Model often feature numerous 
additional interacting heavy particles. The phenomenology of supersymmetric (SUSY) theories, for example, is 
characterized by sparticle production and cascade decays, which lead to many-particle final states and scattering 
amplitudes with complex resonance structure. In order to extract the additional Lagrangian parameters of an 
extended theory from collider data, theoretical predictions are required that match the experimental accuracies. In 
theoretical calculations production and decay stages can be factorized by means of the narrow-width approximation 
(NWA), which effectively results in on-shell intermediate states. Its main advantage is that sub- and nonresonant 
as well as nonfactorizable amplitude contributions can be neglected in a theoretically consistent way, resulting in 
significant calculational simplifications at tree and loop level. For these reasons, the NWA is employed in nearly 
all studies of BSM physics. We note that it is implicitly applied whenever branching ratios are extracted from 
scattering cross sections. A reliable NWA uncertainty determination is therefore crucial. Given the width T and 
mass M of an unstable particle, the uncertainty of the NWA is commonly estimated as 0(T/M). With T/M 
frequently < 2%, its uncertainty is expected to be small in comparison to, for instance, QCD corrections. 

Recently, two circumstances have been observed in which the NWA is not reliable: the first involves decays 
where a daughter mass approaches the parent mass [84], the second involves the convolution of parton distribution 
functions with a resonant hard scattering process [85] . We are thus motivated to investigate when and why the NWA 
is not appropriate in the context of cascade decays in the Minimal Supersymmetric Standard Model (MSSM). We 
first consider a typical example, namely gv,L production at the LHC, i.e. in proton-proton collisions at 14 TeV, 
with the subsequent cascade decay g — > s L s and s L — ► \i c at the SPSla' benchmark point [86, 21] in the MSSM 
parameter space. Phenomenologically, to consider a squark decay into the LSP candidate \i would be more 
natural, but the resulting complete Feynman amplitude features a complicated resonance structure whose study we 
leave to future work. Even for the gluino decay chain considered here, interference arises from g — ► (c* L — ► \i c - 
However, it does not exceed the expected NWA uncertainty and can therefore be neglected. We focus on off-shell 
effects for the resonant sl state (with M = 570 GeV and T = 5.4 GeV at SPSla') and hence treat the chargino 
as stable and the gluino in NWA with spin correlations. As shown in Fig. 1, the NWA error substantially exceeds 
the expectation of T{sl)/M{sl) < 1% when the strange squark mass approaches either the chargino or gluino 
mass of 184 and 607 GeV, respectively. 2 Note that the region where the NWA is inappropriate is not restricted to 
mass configurations where the Breit-Wigner shape is cut off kinematically, i.e. where M(sl) — M(xi) ^ r(sjj 
or M(g)-M(s L )<T(s L ). 

7.2 Resonant 1 — > 3 decays in the MSSM 

The example in Sec. 7.1 suggests resonant 1 — > 3 decays as smallest unit that features the amplified off-shell 
effects. Giving type (4-momentum, mass) for each particle, we define a resonant 1 — > 3 decay by 

T 7 (Pj,M 7 ) ->Ti(pi,mi), T(q,M) and T(q, M) -» T 2 (p 2 , m 2 ), T 3 (p 3 ,m 3 ) . (1) 

The width of the intermediate particle with momentum q is T. Type can be scalar (S), fermion (F) or vector boson 
(V). In the MSSM, 48 generic processes exist and are identified with type codes TiT\T-TT-fr 3 . For each process 
we have systematically scanned the MSSM parameter space for the maximum deviation \R\ of off-shell (r o ff- s heii) 
and NWA (Tnwa) decay rate predictions, where R = (r o ff. s h e ii/rNWA — 1)/(T/M). Note that in R, coupling 
constants typically cancel with the exception of the relative strength of the chiral components of SFF and VFF 
vertices, which has been varied in addition to the masses and width. From this survey [88] we conclude that large 
deviations \R\ do not occur for configurations with a resonance mass that is very far from kinematical bounds. The 
NWA exploits that in the limit T — > the squared propagator D(q 2 ) = [(q 2 — M 2 ) 2 + (MT) 2 ]~ 1 is asymptotically 
equal to 27rif NW A<5(<? 2 - M 2 ) with if NWA = 1/(2 A/T) = D(q 2 ) dq 2 /(2ir). The Breit-Wigner shape is thus 
effectively integrated out. The origin of unexpectedly large deviations for configurations where kinematical bounds 
are outside the resonance region is that the q 2 -dependence of the residual integrand significantly distorts the peak 
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Fig. 1: The accuracy of the NWA cross section normalized to the conventionally expected uncertainty is shown for gu^ 
production at the LHC followed by the cascade decay g — > sls and sl — * Xi c m me MSSM at SPSla' for a variable strange 
squark mass that approaches the chargino mass (left) and the gluino mass (right). 



and tail of D(q 2 ). We find that the effect is most pronounced for the decay process SSS-SSV. We thus use it to 
demonstrate the distortion. With mi = = 0, the q 2 -integrand is given by 

\}~Mj)\ 1 ~~f~)\ m| ) (q 2 - M 2 ) 2 + M 2 r 2 ' (2) 

The 1st- and 2nd-stage decay PS elements contribute the first and second factor, respectively. The 2nd-stage decay 
matrix element gives the third factor. When m 2 < M 2 the second and third factor effect a strong deformation 
of D(q 2 ), which, together with the resulting large deviations, is displayed in Fig. 2. The deviation grows with 
increasing power of the deforming factors. When M approaches the lower kinematic bound, \R\ is sensitive to the 
type of the 2nd-stage decay, which determines the power of the factor that deforms the Breit-Wigner peak. While 
this factor enhances the Breit-Wigner tail, the factor of the lst-stage decay suppresses it. And vice versa for the 
upper bound. We find stronger effects for SSV, VSV, FFV, VVV and SVV than for FSF, SFF, VFF, VSS and SSS 
vertices. Using our generic results we probe resonant 1^3 decays at SPS benchmark points [21]. Decays with 
larger deviations are shown in Table 1 . 

Affected decays generally have a small branching ratio due to similar-mass configurations. For example, 
BR = 1.3% for the decay mode xt ~* Xi u d at SPSla. It proceeds via the intermediate states W + (resonant), 
and d* L (nonresonant). The resonant W + contribution with R-T/M = 13% (see Table 1) induced by the mass ratio 
(mi + M)/Mi — 0.975 dominates, and the off-shell prediction including the nonresonant contributions deviates 
by about 11% from the NWA prediction. Since the lst-decay stage is not affected by QCD corrections, this error 
is particularly significant. For a detailed discussion of effects at SPS points including cascade decay segments we 
refer to Ref. [88]. 



7.3 Conclusions 

When the NWA is applied to decay chains the approximation error will exceed order T/M for mass configurations 
in an extended vicinity of segment kinematical bounds due to a significant distortion of the Breit-Wigner peak 
and tail, which is effected by the g 2 -dependence of the phase space elements and residual matrix elements. In 
phenomenological studies of affected models, fully off-shell tree-level Monte Carlos [87] should thus be used even 
though it requires more computing resources. For decay processes involving strongly interacting particles QCD 
corrections are known to be large [89, 90, 91] and need to be taken into account. For this purpose, a suggestive 
NWA improvement is proposed in Ref. [92]. We have chosen the MSSM to illustrate how large off-shell effects 
can occur in extended models, but emphasize that the effects do not depend on SUSY. 
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Fig. 2: Resonant 1—^3 decay SSS-SSV (see main text) with Y/M = 0.01: The graph displays the g 2 -dependence of the 
Breit-Wigner (dashed) that is integrated out in the NWA and of the complete integrand of Eq. 2 (solid) for 7713 = M — 3T. The 
contour plot shows R, the resulting off-shell-NWA deviation in units of T/M, as function of 7713 and M with mi = mi = 0. 
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0.976 
0.975 
2.49 
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1.11 
1.19 
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Table 1: R, the off-shell-NWA deviation in units of T/M, for resonant 1 — > 3 decays at SPS benchmark points. 
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8. Supersymmetric corrections to Mw and sin 2 l w in mSUGRA 1 
8.1 Introduction 

Specific patterns of supersymmetry (SUSY) breaking provide relationships between various sparticle masses. In 
mSUGRA for instance, the GUT-scale scalar masses are set to m , the gaugino masses Mi/ 2 and the universal 
scalar SUSY breaking trilinear coupling Aq. These degeneracies are broken by renormalisation group effects 
between the GUT scale and Mz- It is by now well known that various regions of mSUGRA parameter space are 
ruled out by direct sparticle search constraints, which place lower bounds upon the sparticle masses. 

Sparticles may appear in loop corrections to electroweak observables, therefore affecting the values of the 
latter as predicted within the Standard Model once a set of independent physical input parameter is chosen, such 
as the electromagnetic couplings, Gf from \i decay, Mz and other SM particle masses. Two such precision 
observables are the W boson mass, Mw, as measured at LEPII and the Tevatron and the effective leptonic mixing 
angle, sin 2 9 l w , derived from the lepton asymmetries measured at LEPI and SLD. The former is related to the 
electric charge e = \/Ama, the weak mixing angle s 2 w = 1 — M w /M^, Gf, and Mz via 

where the parameter Ar is a model dependent quantity which accounts for all higher order corrections to the 
muon decay (this includes self energies, vertex and box corrections in a given model, see ref. [93] for a recent 
discussion in the context of the MSSM). The Mw value which solves the above relation constitutes the model 
specific prediction of Mw for a given fixed set of Standard Model and new physics parameters. Similarly, one can 
express the effective leptonic mixing angle as 

sin 2 ^ = 4(l + AK), (2) 

where the model dependent higher order corrections to the leptonic Z boson decay, Z — ► 11, enter via An (details 
concerning sin 2 9 l w in the general MSSM can be found in ref. [94]). The prefactor s 2 ^ = 1 — M w jM\, and thus 
also sin 2 6 l w , is furthermore sensitive to radiative corrections via Mw- It is convenient to split the MSSM higher 
order corrections into Standard Model and SUSY type contributions [93, 94] 

Ar = Ar SM | M|M=Mh + Ar SUSY , An = An SM \ M%M=Mh + A K SUSY , (3) 

with the Standard Model Higgs boson mass M^ M set to the lightest MSSM Higgs mass M/j. Direct search con- 
straints put lower bounds upon sparticle masses, limiting the size of the SUSY contributions, which are propagator 
suppressed by large sparticle masses. 

Empirical constraints on M w and sin 2 6 l w are very tight: they are taken here to be [95, 96] 

M w = 80.398 ± 0.027 GeV, sin 2 6 l w = 0.23153 ± 0.000175. (4) 

One may ask how large the SUSY contributions to the electroweak observables are, given the current strong 
constraints upon sparticle masses from LEP2 and the Tevatron. If the SUSY contributions are much smaller than 
experimental errors upon the relevant observables, then there is no need to include them in any fit. Here, we 
use results from a previous fit of mSUGRA to dark matter and other indirect data (including the direct search 
constraints) in order to see how big the SUSY contribution to Mw and sin 2 9 l w may be. 



8.2 The fits 

In refs. [95], multi-dimensional fits to mSUGRA were presented using the SOFTSUSY2 .0.10 [97] spectrum 
calculator and the micrOMEGAsl .3.6 [98] dark matter code. m , A , tan/3, M 1 / 2 , m t , mt,, a s (Mz) and a 
were all scanned simultaneously using the Metropolis algorithm in a Markov Chain Monte Carlo technique. It was 
assumed that the WMAP-constrained relic density of dark matter QoAih 2 consisted entirely of the lightest neu- 
tralino, which is stable by the assumption of R-parity. The following data were included in the fit: the anomalous 
magnetic moment of the muon, (g — 2) M , BR(b — > sj), Tevatron BR(B S — > constraints, Q,£>Mh 2 , Mw, 

sin 2 8 l w , LEP2 Higgs constraints as well as other constraints upon sparticle masses from direct searches. Data on 
m t , rrib, a s (M z ) and a were also included in the likelihood. We refer the reader to ref. [95] for the details. We 
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Fig. 1: Distributions to the SUSY contribution to the electroweak observables Mw and sin 2 6 W from mSUGRA fits. The 
posterior probability distribution function (pdf) is shown in red (dark) whereas the profile likelihood is shown in green (light). 
Each histogram has an arbitrary normalisation. Horizontal bands display the one-sided 95% Bayesian credibility (condidence 
level) regions for the posterior pdf (profile likelihood). 



note here that Mw an d sin 2 l w were determined using am embryonic version of the SUSYPOPE [99] code, which 
is a state-of-the-art MSSM calculation [93, 94] of the electroweak observables 2 . SUSYPOPE is also capable of cal- 
culating the Standard Model value prediction for Mw given other Standard Model input parameters and the Higgs 
mass. A list of ~500 000 weighted mSUGRA parameter space points, with weighted frequency proportional to 
their combined likelihoods, was the result of the fit. This set of points is referred to as KISMET (Killer Inference 
in Susy Meteorology). The points are presented for public use on URL 

http: //users . hep forge .org/ "allanach/benchmarks /kismet . html 

In ref. [95], the fits were presented in several ways: the ways relevant to our discussion here will be the 
frequentist fashion (utilising the profile likelihood) and a Bayesian fit with flat priors in the inputs listed above. 
While the profile likelihood takes into account only the best fit parameters, the Bayesian fit includes volume 
effects. Volume effects take into account the volume of the probability distribution: for instance, a region which 
has a very large volume in marginalised (or averaged) parameter directions but a less good fit can still have an 
appreciable effect on the marginalised posterior probability distribution. However, as is well known, the Bayesian 
interpretation is dependent upon the subjective prior choice unless the data are plentiful and precise. As was 
pointed out in ref. [95], the indirect data used in the fits are currently not plentiful and precise enough. We therefore 
display both methods and use the difference between the two as an indication of the size of inherent uncertainty in 
our interpretation of the fits. 

8.3 The W mass and the weak mixing angle 

We take the KI SMET points and re-weight them, taking off the likelihood contributions from Mw and sin 2 6 l w . By 
marginalising against these two variables, we may then examine what size of SUSY contribution to each is expected 
from fits to the other indirect data. In order to calculate the one-sided 95% limits upon the magnitude of the SUSY 
corrections to each observable, we re-bin in terms of | AJ1% | and | A sin 2 8 l w \ . AMw is the difference between the 
mSUGRA prediction and the Standard Model one with a s , a, m t , mi, and Mh identical to those associated with 

2 The SOFTSUSY determination of Mw and sin 2 9 l m is not to a sufficient accuracy, given the tiny empirical errors upon their 
measured values. 
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the mSUGRA point in question. The posterior pdf with flat priors gives us a 95% Bayesian credibility interval of 

\AM W \ < 19MeV, |A sin 2 6 l J < 12 x 10~ 5 . (5) 

The 95% upper confidence limit in the frequentist interpretation is obtained in each case from the profile likelihood: 

\AM W \ < 70MeV, |A sin 2 9 l J < 30 x 1(T 5 . (6) 

Thus, the frequentist constraints are somewhat (two and a half times) more relaxed than the Bayesian constraints 
with flat priors. These numbers are to be compared with the empirical uncertainties quoted in Eq. 4 of 17.5 MeV 
and 27 x 10 -5 , respectively. Even the tighter 95% Bayesian constraints are around the same values and so we 
conclude that the SUSY contribution to the likelihood cannot be neglected. 
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9. LHC and the muon's anomalous magnetic moment 

9.1 Introduction 

The strongest hint for a TeV-scale modification of the Standard Model originates from the anomalous magnetic 
moment of the muon. This parameter has been evaluated both in experiment and in theory to unprecedented 
precision. We use the experimental value [100] 

a|f xp) = {g - 2)^/2 = 116592080(63) x 10~ n . (1) 

In contrast, the Standard Model prediction is smaller [101, 102] 

a (SM) = n 6 59i7 8 5( 6 i) x iq- 11 , (2) 

where we use the number using e+e~ data, as the discrepancy with tau data is still large [103]. This difference 
corresponds to a 3.4er deviation between theory and experiment, suggesting a new-physics effect. Supersymmetry 
provides a particularly attractive explanation of this discrepancy. Contributions from the supersymmetric partners 
of the muon, the muon-neutrino and the gauge and Higgs bosons modify the Standard-Model prediction. The 
masses of the new particles responsible for this signal should be of the order of several hundred GeV, a mass range 
well accessible to the LHC [104, 105]. Therefore, LHC will be able to test the hypothesis that this discrepancy is 
caused by TeV-scale supersymmetry. 

However, the benefit of such a test will go in both ways: after the LHC will have measured observables 
like the masses of supersymmetric particles or kinematic edges involving such particles [106, 61], the question will 
arise what the fundamental parameters of the Lagrangian are [3 1 , 32] . Since (g — 2) at leading order is proportional 
to tan (3 it includes useful information [107], which can significantly improve the extraction of tan j3. 

As an example we use the experimentally well-studied parameter point SPSla [21]. Its theory prediction 
for a M is a^ SPSla - ) = a|f M ^ + 282 • 10 -11 . This leads to a deviation from the experimentally observed value 
of Aa M = —13 • 10 -11 , which is well below the experimental error bounds. Therefore, we can safely use the 
experimental value without further modifications. 

Our analysis uses the parameter extraction tool SFitter [31] where we have added the necessary modules 
to calculate the anomalous magnetic moment. To obtain a good balance between precision and required time to 
perform the scans we use the one-loop expression for (g — 2) with additional leading two-loop QED-logarithms. 
Note that in this study we are not mainly interested in the best-fitting MSSM parameter point, but in the errors of 
the MSSM parameters, so this simplification is appropriate. 

9.2 Weak-scale MSSM analysis 

The determination of the weak-scale MSSM Lagrangian at the LHC is clearly preferable to tests of SUS Y-breaking 
assumptions, as long as we have enough information available at the LHC. The extracted model parameter can then 
be run to a higher scale, to test for example unification patterns [31]. However, some model parameters can be 
fixed, because neither LHC nor (g — 2) will include any information on them. Properly including the top-quark 
mass as a free parameter we assume a 19-dimensional weak-scale MSSM parameter space listed in Table 1 . As we 
have shown in Ref. [31], even this reduced parameter space cannot be determined completely at the LHC. In the 
parameter point SPSla, we for example find an eightfold degeneracy in the gaugino-higgsino sub-sector. Because 
only three of the neutralinos and none of the charginos can be observed at the LHC, the connection between their 
masses and the model parameters Mi, M 2 and [i is not unique. In addition the sign of \i is not determined by LHC 
data alone [108]. 

This is where the anomalous magnetic moment of the muon adds important information. First, the deviation 
from the Standard Model prediction is proportional to the sign of \i, the parameter which couples the two Higgs 
superfields in the superpotential. Including (g — 2) data will clearly favor one sign of \i, namely the [i > 0, thereby 
reducing the degeneracy by a factor of two. 

When reconstructing the fundamental parameters of the Lagrangian the central values have to be accom- 
panied with the correct error bars. There are three different types of experimental errors on the observables: a 
statistical error and the two (correlated) systematic errors for the jet and lepton energy scales. All experimental 
errors are Gaussian shaped. In addition, we include flat theory errors of 1% for all colored particle masses and 
0.5% for all others. For (g — 2) we use the values given in Eqs.(l, 2). The convolution of these errors is described 
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only experimental errors 


including flat theory errors 


SPSla 




LHC 


LHC(g)(g-2) 


LHC 


LHC<g>(g-2) 




tan (3 


9.8± 2.3 


9.7± 2.0 


10.0± 4.5 


10.3± 2.0 


10.0 


Mi 


101.5± 4.6 


101. 1± 3.6 


102.1± 7.8 


102.7± 5.9 


103.1 


M 2 


191.7± 4.8 


191.4± 3.5 


193. 3± 7.8 


193.2± 5.8 


192.9 


M 3 


575.7± 7.7 


575.4± 7.3 


577.2± 14.5 


578.2± 12.1 


577.9 


M fL 


196.2±O(10 2 ) 
136.2±O(10 2 ) 


263.4±O(10 2 ) 
156.8±O(10 2 ) 


227.8±O(10 3 ) 


253.7±O(10 2 ) 
134.1±O(10 2 ) 


193.6 


Mf R 


164.1±O(10 3 ) 


133.4 


M H 


192.6± 5.3 


192.3± 4.5 


193.2± 8.8 


194.0± 6.8 


194.4 


Ma R 


134.0± 4.8 


133.6± 3.9 


135.0± 8.3 


135.6± 6.3 


135.8 


Me L 


192.7± 5.3 


192.2± 4.5 


193. 3± 8.8 


194.0± 6.7 


194.4 


Me R 


134.0± 4.8 


133.6± 3.9 


135.0± 8.3 


135.6± 6.3 


135.8 




478.2± 9.4 


476. 1± 7.5 


481.4± 22.0 


485.6± 22.4 


480.8 




429.5±0(1O^) 


704.0±O(10 ) 


4l5.8±O(l0^j 


439.0±O(10" : ) 


408.3 


M b R 


501.2± 10.0 


502.4± 7.8 


50l.7± 17.9 


499.2± 19.3 


502.9 


m Sl 


523. 6± 8.4 


523. 0± 7.5 


C r\ A f I 1 A C 

524. 6± 14. 5 


525.5± 10.6 


526.6 


M q R 


506.2± 11.7 


505.8± 11.4 


507.3± 17.5 


507.6± 15.8 


C AO 1 

508.1 


A T 


fixed 


fixed 


fixed 


fixed 


-249.4 


A t 


-500.6± 58.4 


-519.8± 64.3 


-509. 1± 86.7 


-530.6± 116.6 


-490.9 


A b 


fixed 


fixed 


fixed 


fixed 


-763.4 


m A 


446.1±O(10 3 ) 


473.9±O(10 2 ) 


406.3±O(10 3 ) 


411.1±O(10 2 ) 


394.9 




350.9± 7.3 


350.2± 6.5 


350.5± 14.5 


352.5± 10.8 


353.7 


m t 


171.4± 1.0 


171.4± 1.0 


171.4± 1.0 


171.4± 0.90 


171.4 



Table 1: Result for the general MSSM parameter determination at the LHC in SPSla. The left part neglects all theory errors, 
the right one assumes flat theory errors. In all cases a set of 20 kinematic endpoints and the top-quark and lightest Higgs-mass 
measurements have been used. In the third and fifth column we include the current measurement of (g — 2). All masses are 
given in GeV. 

in Ref. [31]. To determine the errors on the model parameter we randomly smear the nominal values for SPSla. 
The corresponding random numbers obey a distribution according to the associated errors. Then we minimize x 2 
for each pseudo-measurement and repeat this procedure 10000 times. The emerging distribution of the parameters 
is simply the result of the correct error propagation. Using a Gaussian fit we then extract the central value and the 
lcr standard deviation of each parameter. 

Table 1 shows the result of our SPSla analysis. For comparison and to make the effect of the additional 
(g — 2) data easily visible, we include the result without (g — 2) data from Tables VIII and IX of Ref. [31]. We give 
results with experimental errors only (columns 2 and 3) and including theory errors (columns 4 and 5). The effect 
of the additional information on the accuracy of the parameter determination is clearly visible. It is particularly 
significant for tan (3, which is not well determined by the measurements of kinematic endpoints at the LHC. The 
best source of information on tan j3 is the light MSSM Higgs mass [109], but this observable strongly relies on the 
assumed minimal structure of the Higgs sector, on the knowledge of many other MSSM parameters, and on the 
estimate of the theory errors due to higher orders. Because of a lack of complementary measurements (for example 
A t ) a change in tan j3 can always be compensated by an appropriate change in other MSSM parameters, leaving 
the value of all LHC observables unchanged. Additional sources of a tan (3 measurement are the production rate 
for heavy Higgs bosons [110] and rare decays like B s — > ^ + ^, which we study elsewhere in this volume, but 
both of them only work for large enough values of tan f3. 

The (g — 2) prediction has a leading linear dependence on tan (3. Therefore, the improvement of the 
tan (3 errors by more than a factor of two can be easily understood. This improved accuracy of tan (3 influences 
those parameters which must be re-rotated when tan/3 is changed to reproduce the same physical observables. 
Correlations and loop corrections propagate the improvement over almost the complete parameter space. 
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9.3 SLHAio 

The link between the main SFitter modules and the (g — 2) module is provided by SLHAio. SLHAio is a library 
which allows for a smooth communication between different programs according to the SUSY Les Houches Ac- 
cord (SLHA) [72] and its extension SLHA2 [111]. With such an interface, for example a (g — 2) code can easily 
be used in a large-scale program like SFitter. 

The separation of actual calculations and SLHA interface creates a simple structure. SLHAio itself has an 
easy instruction set with a concept similar to the XML Path Language (XPath) [112]. Each data field in SLHA is 
associated with a leaf in a tree. The leafs are accessed and manipulated with a path. The path itself is similar to the 
path used for navigation in file systems. 

The tree has no predefined structure. Its size is determined by the amount of information stored. It can 
of course grow beyond the fields defined in SLHA and SLHA2. This feature we use for the smuon mass matrix, 
where SLHA restricts itself to third-generation particles. There are two data types for each field: string and double. 
Conversions are done automatically. Setting a leaf with a string and afterwards reading a double is possible, as long 
as the string can be converted. The library assures the highest precision possible. Strings which are never changed 
through doubles remain strings, because this representation has the highest precision. Conversions to strings are 
done in the format defined in SLHA. One part of SLHAio is a templated container class for matrices. This class is 
fully integrated into SLHAio, so matrices can be read, stored and printed with a single command. 

As discussed above, the (g — 2) code can be used stand-alone and within SFitter. When used with SFitter, 
data is shared directly via SLHAio. This is a huge increase in performance, because no files need to be written and 
doubles do not need to be converted to strings and back again. It also means that once the tree is set up, parameter 
changes are as fast as the access to a double pointer. At this stage, no SLHAio functions are involved. While 
SLHAio is used in the SFitter extraction tool, it will in the future become publicly available. 

9.4 Conclusions 

Supersymmetry provides a particularly convincing explanation for the currently observed 3.4er discrepancy be- 
tween the experimental value for the anomalous magnetic moment of the muon and its Standard-Model prediction. 
If this signal proves to be correct, there exist new particles in a mass range accessible to the LHC. When we attempt 
to reconstruct weak-scale MSSM parameters from LHC observables, the numerical value of (g — 2) provides an 
attractive additional handle on the MSSM parameters. Using the parameter point SPSla as an example, we have 
shown that (g — 2) essentially determines the sign of /i in the weak-scale MSSM, cutting the number of discrete 
solutions obtained in the general MSSM in half. In addition, many of the parameter errors are reduced, most 
notably the error on the extraction of the notoriously difficult parameter tan 0. 
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10. Towards combining cascade decays and flavor physics 
10.1 Determining supersymmetric parameters 

For many years it has been known that the kinematics of cascade decays is particularly well suited to extract 
the masses of the particles involved [106, 61]. More recently, we have seen how these masses can be used to 
determine the (running) TeV-scale Lagrangian parameters, with the ultimate goal of evolving these parameters to 
higher energy scales and extracting information on the mechanism of supersymmetry breaking. SFitter [31] and 
Fittino [32] are two computer tools specifically designed to determine TeV-scale supersymmetric parameters with 
the proper experimental and theory errors. 

Most studies of cascade decays are based on the decays of (five) light-flavor squarks and gluinos. Hence, 
for example in SPSla we would have convincing control over these squark masses and over the neutral gaugino 
masses. A crucial and yet notoriously hard parameter to extract is tan (3, both in the neutralino/chargino sector 
and in the Higgs sector. In most analyses, we rely on the light MSSM Higgs mass for information on tan (3 [109]. 
This extraction depends on a large number of supersymmetric parameters, on the strict MSSM assumption, and on 
a reliable estimate of the theory errors. Elsewhere in this volume we show how a measurement of (g — 2)^ can 
be used to determine tan (3 from the lepton sector. In a former Les Houches project it was shown that tan (3 can 
be extracted from the production rate of heavy MSSM Higgs bosons, which for large values of tan (3 is typically 
proportional to to 2 tan 2 (3. Combining all errors entering the cross-section measurements this study predicts a 
total error of 12% to 16% on the bbA/bbH Yukawa coupling, which is proportional to tan f3 [110]. 

Another strategy for an indirect tan/3 measurement are flavor-physics observables. For example the rare 
decay rate for B s — > 11 is proportional to tan 6 (3/m\. This steep behavior makes it a prime suspect to extract 
tan (3 [113]. A major problem of such an extraction is the correct estimate of the theory error on the observable. 
The second problem is the dependence of the effective bs{h, H, A} couplings on the stop and chargino masses 
appearing in the loop. In this article we briefly report on a pre-study done for SFitter, to give a first estimate if 
these two problems will leave channels like B s — > 11 promising candidates to be included in the SFitter set of 
observables. 

Unfortunately, the usual parameter point SPSla with tan /3 = 10 is not well suited to study the determina- 
tion of tan (3. Even if there should be a sensitivity from a measurement of B s — > 11, it is unclear if we will observe 
any of the heavy Higgs bosons at the LHC. We therefore modify this parameter point in the direction of SPSlb, 
simply choosing a range of larger tan (3 values. For simplicity we assume that the set of cascade observables is not 
altered by this change, including the sbottom mass determination. While this assumption might be quantitatively 
naive, it will serve our purpose of estimating the odds of combining different sources of information on tan (3. 



10.2 Combining flavor and cascades 

FCNC processes involving down-type quarks in the Standard Model are both highly suppressed and sensitive to 
the mass and couplings of the top. This is because, beyond their loop suppression, the unitarity and hierarchical 
structure of the CKM matrix and the hierarchy m t . w ^ fn UjC , entail strong GIM cancellations between the light- 
flavor loop contributions. Turning around this argument, they are sensitive probes of new-physics effects. In 
particular when the supersymmetric flavor structure is (close to) minimally flavor-violating [1 14, 1 15, 1 16], they 
provide a handle on stop and chargino masses. Of special interest are FCNC mediated by neutral Higgs exchange, 
which exhibit a double enhancement: first, they involve the large bottom Yukawa coupling y^ oc yf M tan j3. 
Secondly, the loop-induced contribution of v u = (H u ) >• Vd to the down-type fermion mass matrix destroys the 
alignment of the mass matrix with the Yukawa couplings and renders the latter flavor-nondiagonal [117, 118, 119, 
120, 121, 122, 123], leading to an additional factor tan (3. For minimal flavor violation, the corresponding FCNC 
Higgs couplings have the form [123, 124, 125] 

2 

£cff D b R a L (t>i Xi V t * b V ts j- -, ■ — % — \T o\n I m ' (!) 

cosp (1 + (e + y\ fey) tan/j)(l + e tan/3) 

where xh = — sin(a — (3) and xa — i- The FCNC couplings of h are suppressed by a factor Xh = cos(a — (3), 
offsetting the tan (3 enhancement. The parameters eo an d e Y parameterize loop-induced "wrong-Higgs" contribu- 
tions to the down-quark mass matrix. The sensitivity to MSSM parameters becomes most transparent in the limit 
v <C Msusy, when 

xlnx y In y 



1 A t 

(Y = 



16 7T 2 y t [i L(l - x) [x - y) (l-y)(y-x) 
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with x = m-^/fj 2 and y — m^/^ 2 . The observables most sensitive to these couplings are B s ^ — > , where 
the tree-level iJ, A exchange contributes at the amplitude level as [123, 124, 125, 126, 127] 

A(B q -» ^/")^ oc -J£« ^ 1 Wt 2 an3/3 . (3) 
cos p m m ^ 

Among the modes accessible at LHCb, ATLAS, and CMS, B s — > /i + yU~ has the largest branching fraction [113], 
which can be dominated by the neutral-Higgs contribution. 

The usual cascade measurements we assume with this simple study include g — > tit — ► t?Xi> 9 ~~ * £ii — > 
*&Xi~> <Zl -» X§9 -» r fl rZg -> XoTTqZ, U -» -> X^9 -> qW\iZq -> qq'q"xiZ, where 

the xj cascade is strictly speaking not necessary for our analysis. In our simple parameter points M = 150, 

= 250, A = -100, (i > 0, and tan/3 = (30, 40) [128], these cascades should be visible. In our Minuit 
fit we include 10 observables listed in Table 1, including BR(B — > /^u) and the edge measurement m t b [129]. In 
addition, we need to include some very basic information on the chargino-stop sector. The measurement of the 
three neutralino masses gives us information on the chargino mass parameters M2 and \i. The left-handed stop 
mass is linked to the left-handed sbottom mass via SU(2). However, the right-handed stop mass as well as the off- 
diagonal entry into the mass matrix are not determined by cascade decays. The latter is dominated by the trilinear 
coupling A t , which enters for example the calculation of the light Higgs mass [130]. Extracting A t , however, 
requires a measurement of the dominant heavy Higgs mass parameter, which again limits us to reasonably large 
values of tan (3. We assume vha to be known either directly or via the charged Higgs mass, similarly to Ref. [110]. 

The modified SPSla parameter point for two example values of tan/3 is specified in Table 2, together with 
the best-fit values and the errors from our Minuit fit. For B s — > fi^i, LHCb alone expects about 100 events at the 
Standard-Model rate after 5 years of running [113]. For our study we assume an integrated luminosity of 10 fb _1 
for the B s sample. The Higgs-mediated contributions always increases the corresponding events number. As a 
consequence, the theory error, which at present ranges around 30%, will soon dominate the total uncertainty, unless 
it can be reduced. It is mainly due to the decay constant fg e , which can be calculated using numerical lattice-QCD 
methods (see Ref. [131] for a recent review). For our study we simply assume a reduction of the error on fs s to 
7%, about half its present value. Such a reduction is commonly believed to be realistic over the next five years. 

For our study, we perform two sets of fits, one ignoring the theory error and one combining it in quadrature 
with the statistical error. A more refined treatment of the theory error is in progress and will use the proper Rfit 
ansatz [43], as implemented in Sfitter [31]. We indeed see that without taking into account the theory error, tan/3 
will be determined to 10% from the combined toy data sample. Including a realistic theory error increases this 
number to 15 • • • 20%. The errors on the remaining parameters, shown in Tab. 2 remains largely unchanged. Slight 
shifts in either direction are at his stage well within the uncertainly on the determination of the error bars, and 
the central fit value for the top-mass parameter seems to be consistently lower than the input value (by roughly 
half a standard deviation). Comparing our error estimates on the mass parameters for example with the SFitter 
analysis [31], we expect the situation to improve for all model parameters once we include a more extensive set of 
measurements and properly correlated errors. 



tan/3 


30 


40 




value error 


value error 


m h 

m t 
m H ± 
m x o 

X2 

TO ± 

Al 

Wig 

BR(B S -» w) 
ma 


112.6 4.0 
174.5 2.0 
354.2 10.0 
98.4 4.8 

183.1 4.7 
353.0 5.1 
182.8 50.0 

607.7 8.0 
7.3 -lO" 9 VN®15% 

404.2 5.0 


112.6 4.0 
174.5 2.0 
307.2 10.0 
98.7 4.8 

183.5 4.7 

350.7 5.1 

183.1 50.0 

607.6 8.0 
3.2-10~ 8 VN®15% 

404.2 5.0 



Table 1: Set of toy measurements. The simple combined (absolute) errors are SPS la-inspired. 
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no 


theory error 


ADD /T?R 


— 1O/0 


no theory error 


ART? /P.T3 


— ±0/0 




true 


best 


error 


best 


error 


true 


best 


error 


best 


error 


tan/3 


30 


29.5 


7 A 

3.4 


29.5 


6.5 


A A 

40 


7n o 
39.2 


A A 

4.4 


39.2 


5.8 


n/r 

M A 


7/m 7 


344.4 


33.8 




5Y.L 


295.5 


304.4 


35.4 




77 Q 
jj.y 


Mi 


101.7 


100.9 


16.3 


100.9 


16.4 


101.9 


101.0 


16.3 


101.0 


16.3 


M 2 


192.0 


200.3 


18.9 


200.3 


18.8 


192.3 


200.3 


20.0 


200.7 


18.9 




345.8 


325.6 


20.6 


325.6 


20.6 


343.5 


322.9 


20.7 


323.3 


20.6 


M 3 


586.4 


575.8 


28.8 


575.8 


28.7 


586.9 


576.0 


28.7 


575.8 


29.0 


M Q L 


494.4 


494.4 


78.1 


494.3 


78.0 


487.1 


487.6 


79.4 


487.5 


78.9 


M tn 


430.0 


400.4 


79.5 


399.8 


79.5 


431.5 


399.2 


86.7 


399.1 


82.6 



Table 2: The modified SPSla point and the errors from the parameter fit for the two values of tan (3 = 30, 40. Dimensionful 
quantities are in units of GeV. For the measurement of BR(B S — * it) we assume either no theory error or an expected 
improvement to 15%, as compared to the current status. 

10.3 Caveats 

Note that the detailed results of this study should not be used at face value. First of all, it is not clear if the stop- 
mass measurement can be achieved in the SPSla parameters point or with an increased value of tan/3. Secondly, 
for the charged Higgs mass we only use a toy measurement. And last but not least, we do not (yet) take into account 
error correlations at this stage. None of these omissions we expect to move the result of a complete analysis into a 
definite direction, but there is certainly room for the final error bars to move. 

This study shows, however, that the parameter tan (3 can indeed be extracted from a combined cascade 
and flavor data sample. Already at this stage we can conclude that the combination of cascade-decay and flavor 
observables will crucially depend on the quality of the theory predictions in the flavor sector. In particular, an 
improved understanding of non-perturbative QCD effects in B s — > U decays is needed to meaningfully exploit 
this highly promising link. From the high-py point of view we also generally see that to measure tan (3 we need 
to improve the analysis of the stop-chargino sector in the classic decay-kinematics analyses. 

The problem with the measurement of tan/3 from the light MSSM Higgs mass or from (g — 2) M or from 
rare B decays is that each of these indirect measurement rely on assumptions about the flavor and Higgs sectors. 
Moreover, these different measurements point to different parameters, not only from a renormalization point of 
view, but also because of large QCD effects distinguishing between them. The more direct extraction from cross 
sections times branching ratios of heavy Higgs bosons is at the same time plagued by large theory uncertainties, due 
to QCD corrections and uncertainties in the bottom-parton picture [132, 133]. If we should indeed find evidence 
for the MSSM in the LHC era we obviously expect a serious jigsaw approach to the tan (3 determination. 

Acknowledgments 

We are grateful to Mihoko Nojiri and Giacomo Polessello for their ongoing encouragement. And of course we 
thank the organizers of the great Les Houches 2007 workshop Physics at TeV-Scale Colliders] 



46 



11. BBN lithium problem consequences at LHC 1 

11.1 Lithium problem 

Recent measurements of the fluctuations of the microwave background radiation allowed for determination of 
the contribution of baryons to critical density f^/i 2 w 0.0224 [36] which is the only free parameter of the stan- 
dard BBN. SBBN previsions are in agreement with experimental estimates of abundances for deuterium and and 
helium 4 but have problem in explaining lithium 7 proliferation, which is about 3 times to high. 

There is no agreement whether stated above discrepancy is a real problem because apparent primordial 7 Li 
abundance is derived only from the observation of so called Spite plateau [134] in low metallicity POP II stars. It 
is believed, that gas present in atmospheres of these very old starts have not changed composition since the BBN 
era. 7 Li is, however, fragile so in principle its depletion could be explained by some stellar evolution model but no 
fully satisfactory model have been proposed [135, 136]. Recently preliminary observation of similar plateau was 
done also for 6 Li [135] which was produced during SBBN but below detectable level. If confirmed, the plateau 
could suggest pre-galactic 6 Li origin corresponding to the primordial abundance at least an order of magnitude 
higher than SBBN predictions. Moreover, since 6 Li is by far more fragile than 7 Li, any model of its destruction 
will aggravate 6 Li problem. 

Although above sketched lithium problems could have standard explanation it is very interesting to note 
beyond standard model solutions that alters BBN. All of them postulate long lived massive particles with lifetimes 
around Iks or more. Their decays allow for late 6 Li production or late 7 Li destruction or both without altering 
abundances of other isotopes. If these particles could be negatively charged than also bound states must be taken 
into account. In the present letter there is no place for review all proposed solutions. We would concentrate on two, 
which are in agreement with cosmological constraints, and could have very interesting consequences for the LHC 
phenomenology. In both, an existence of long-lived stau is postulated. In the first [137] stau lifetime is of the order 
of Iks. The solution for both lithium problems is found within CMSSM for stau mass of the order of 1 TeV which 
is, unfortunately, out of reach at LHC. However, solution of only 6 Li problem is possible for stau mass around few 
hundreds GeV. The second solution of both lithium problems is possible for very long stau lifetime of 1 Ms and 
stau mass of about 300 GeV [136]. 

11.2 Possible discovery at LHC 

Both ATLAS and CMS experiments have developed strategies to look for charged massive particles (CHAMP) if 
they decay lengths exceed detector sizes. 

Methods tested on full detector simulation are based on TOF measurements in the muon systems (CMS 
drift tubes, ATLAS drift tubes and RPCs) or specific ionization measurement in the tracker (CMS). Since in 
both experiments at least two independent measurements are performed it is possible to evaluate misidentification 
probabilities directly from data. If simulated performances will be confirmed almost background free selections 
with efficiency of the order of 10% could be designed. After collecting 10/fb of data this allows for the discovery 
if the cross section exceeds lOpb which corresponds to stau mass around 300 GeV [138]. 

This provoke natural question about lifetime measurement of such CHAMP. Any estimate of it could be of 
crucial importance. The problem is, however, that interesting range of lifetimes is above 100 s. It is obvious that 
decays in flight are not only by far inefficient but also insensitive. 

However, the first proposition in this direction was made already a decade ago [139]. The point was about 
using CMS muon system as a late electromagnetic calorimeter (so called /iCAL) in the following way. If significant 
fraction of energy release is electromagnetic and if the decay happens inside iron yoke at a distance to the next muon 
station equivalent to few radiation lengths, then developing electromagnetic cascade causes large accumulation of 
hits in the station. By design the cascade ends in the next yoke section. So the signal is large accumulation of hits 
in one muon station. Despite the fact, that this proposal was made in the context of detection of decaying in flight 
neutralinos, it could be used also for decaying staus, but to measure long lifetimes these staus must stop inside the 
yoke. Unfortunately only small fraction of staus will do that. Larger fraction will be stopped in the concrete and 
rocks of the cavern. 

There were proposals to drill out the cavern walls to recuperate the part with stopped CHAMP or to install 
water tanks for CHAMPs capture, but if CHAMPs are staus than there is also much simpler solution. 17% of 
stau decays produce muons. These muons could be detected not only for staus stopped inside detector but also 
for staus stooped in the cavern walls if muon is released backward. The first possibility is rather hopeless for the 
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ATLAS detector because of air muon system but the second could be more efficient for ATLAS than for CMS 
because ATLAS is bigger and is closer to the cavern wall and so, it will have larger angular size when seen from 
the CHAMP decay point. 

Although it is straight forward to estimate for a given model a probability for a muon from the stopped stau 
decay to cross again the detector full simulation is needed to know detector answer to such muon. In principle 
these muons are very similar to cosmic muons. The only difference is a homogeneous distribution of incoming 
directions in contrast to top bottom directions expected for cosmics. 

Without full detector simulation (which already started, but no official results of it have been released) it 
is difficult to know if the sensitivity could be sufficient for measuring long lifetimes of stau. However, no such 
measurement is possible without cosmic trigger. Another problem is tracking of low momentum staus, which will 
go outside 25 ns window. Although triggering on such particles is not possible its offline recovery is not hopeless 
because drift tubes, thanks to their operation mode, remember data form many bunch crossings. 

It is important to note, that cosmic trigger, even if possibly harmful for normal LHC operation, should be 
studied in detail, to be switched on if long lived CHAMPs will discovered at LHC. 

11.3 Conclusions 

It was underlined that lithium BBN problem could be solved by long lived CHAMPs, next, that such particles, if 
not to heavy, could be discovered at the LHC and that the most efficient way of measuring their very long lifetime 
is to design cosmic trigger for LHC detectors. 
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12. Precision measurements of the stop mass at the ILC 1 
12.1 Introduction 

Supersymmetric particles are likely to be produced and observed in high-energy proton-proton collisions at the 
LHC. However, it will be difficult to confirm their identity as superpartners of the known Standard Model particles 
and to measure their properties precisely. For this, one needs experiments at a linear c+c~ collider such as the 
proposed ILC at ^/s = 500 GeV. The importance of scalar top studies has been emphazised in the '2005 Les 
Houches' proceedings [140]. This work extends these studies. 

An experiment at the ILC will be able to make many precise measurements from which particle properties, 
and ultimately, the outlines of a particle physics model may be inferred. Due to the high statistical precision 
expected at the ILC, the optimization of the systematic errors is of particular importance. We have studied one 
specific example, the extraction of the mass of a scalar top quark from cross-section measurements near threshold. 
We have devised a method which reduces most systematic uncertainties and leads to a potentially very accurate 
measurement of the stop quark mass. This method is general and could be applied to other particles that are 
pair-produced in an e + c~ collider. 

The method relies on the comparison of production rates at two different center-of-mass energies, and 
knowledge of how the cross-section varies as a function of y/s and the mass of the particle. 

We have chosen the case of a light scalar top with a mass not much higher than the mass of the lightest 
neutralino since production of this particle was already extensively studied in an ILC context. It was concluded 
that a conventional approach to the measurement of the stop quark mass culminated in an uncertainty of about 
1 GeV [141]. The new method improves substantially on this result. The presented results are preliminary and 
being finalized [142]. 

For this analysis, we have performed realistic simulations of the signal and backgrounds, and used two 
techniques to separate the signal from the background. The first technique is based on conventional sequential 
cuts, while the second employs an Iterative Discriminant Analysis (IDA). Furthermore, the hadronization followed 
by fragmentation of the stop has been included and we have carefully studied the systematic uncertainties arising 
from this and other sources. 

There are theoretical motivations for studying a light stop quark with a small mass difference. Specifically, 
we evoke a scenario within the Minimal Supersymmetric extension of the Standard Model (MSSM) which is able 
to explain the dark matter density of the universe as well as the baryon asymmetry through the mechanism of 
electroweak baryogenesis [141]. 

A small mass difference between the stop and the lightest neutralino can help to bring the dark matter relic 
density into the observed region [36, 143] due to co-annihilation between the stop and the neutralino. For this 
mechanism to be effective, the typical mass difference is rather small, Am — m~ ti — m^o < 30 GeV [144]. The 

dominant decay mode of the stop is ti — > cx?> resulting in a final state with two soft charm jets and missing 
energy. 

Previous methods to determine the scalar top 
quark mass were discussed for the SPS-5 benchmark 
(m^ = 220.7 GeV) [145] and results are summarized 
in Table 1 . For the cosmology motivated benchmark 
with mr = 122.5 GeV and m^o = 107.2 GeV, an 



Method 


Am;, (GeV) 


luminosity 


Polarization 


0.57 


2 x 500 fb" 1 


Threshold scan 


1.2 


300 ftr 1 


End point 


1.7 


500 ftr 1 


Minimum mass 


1.5 


500 ftr 1 



experimental precision of Am; = ±1.0 GeV was ob- 

. . , ,,,,, , , . or u' i a- »u Table 1 : Comparison of precision for scalar top mass determi- 
tained [141], and about ±1.2 GeV including theoreti- ft- r 

cal uncertainties. The following study investigates the nation for the SPS " 5 benchmark (m h = 220.7 GeV). 
same signal scenario and it is based on the same back- 
ground reactions and event preselection. 



12.2 Mass determination method 

This method proposes to derive the stop mass from measurements at two center-of-mass energies, one measuring 
the stop production cross-section near the threshold (th), and the other measuring it at a center-of-mass energy 
where the cross-section has approximately a peak (pk). Using both measurements leads to a cancellation of sys- 
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tematic uncertainties in the mass determination. A parameter Y is defined as 

N th - B th 



Y = 



N pk - B 



P k 



fth £th Ah 
Opk £pk £pk 



(1) 



where N is the total number of expected events after event selection and B the number of corresponding back- 
ground events, a is the stop production cross-section, e the selection efficiency, and C the luminosity. The center- 



of-mass energies 260 and 500 GeV have been chosen, 
sensitive to the stop mass. 

In this study we assume that the ILC will op- 
erate primarily at y/s = 500 GeV with a total lumi- 
nosity of C = 500 fb -1 , and a small luminosity of 
C = 50 ftr 1 will be collected at y/s = 260 GeV. 
Table 2 summarizes the expected production cross- 
sections. The detector response was modeled with the 
SIMDET package [146] including the LCFI [147] ver- 
tex detector concept. 

The relation of the observable Y and the stop 
mass is used to determine the stop mass with preci- 
sion. For example a variation of Y by 3% in a realistic 
scenario would lead to an uncertainty of the stop mass 
Ami = O-^ GeV as illustrated in Fig. 1. 

12.3 Sequential-cut analysis 

In order to cancel the systematic uncertainties to a large 
extent with the described method, the same sequential 
cuts are applied for the y/s — 260 and 500 GeV analy- 
ses. Details of the event selection are given in Table 3 
and the results are given in Table 4. 



Near the threshold, the production cross-section is very 
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0.13 




0.12 - 



0.11 - 



0.1 - 



0.09 

121.5 122 122.5 123 123.5 

stop mass (GeV) 
Fig. 1: Example of mass uncertainty derivation from the uncer- 
tainty of the observable Y. 



Table 2: Cross-sections for the stop 
signal [148] and Standard Model 
background processes for y/s = 
260 GeV and y/1 = 500 GeV and dif- 
ferent polarization combinations. The 
signal is given for a right-chiral stop 
of m t - = 122.5 GeV. Negative polar- 
ization values refer to left-handed po- 
larization and positive values to right- 
handed polarization. 



Process 


a (pb) at y/s = 


260 GeV 


a (pb) at y/s = 


500 GeV 


P(e-)/P(e + 


) 0/0 


-.8/+.6 


+.8A.6 


0/0 


-.8/+.6 


+.8/-.6 


tit* 


0.032 


0.017 


0.077 


0.118 


0.072 


0.276 


W + W~ 


16.9 


48.6 


1.77 


8.6 


24.5 


0.77 


zz 


1.12 


2.28 


0.99 


0.49 


1.02 


0.44 


Wev 


1.73 


3.04 


0.50 


6.14 


10.6 


1.82 


eeZ 


5.1 


6.0 


4.3 


7.5 


8.5 


6.2 


qq,q/t 


49.5 


92.7 


53.1 


13.1 


25.4 


14.9 


tt 


0.0 


0.0 


0.0 


0.55 


1.13 


0.50 


2-photon 

Pt > 5 GeV 


786 






936 







Variable 



Table 3: Selection cuts for y/s = 260 
and 500 GeV. Also listed are the selec- 
tion efficiencies optimized for right- 
chiral stop quarks. 



y/s = 260 GeV 



y/s = 500 GeV 



number of charged tracks 
visible energy E v i s 
event long, momentum 
event transv. momentum pr 
thrust T 

Number of jets Nj cts 
extra-jet veto 

charm tagging likelihood P c 
di-jet invariant mass rrijj 



5 < V tracks < 25 
0.1 < E vis /y/s < 0.3 
|pi,/ptot| < 0.85 
15 < Pt < 45 GeV 
0.77 < T < 0.97 

N icts > 2 
E jct < 25 GeV 
P c > 0.6 
ri| < 5500 GeV 2 or 
> 8000 GeV 2 



m j 2 j 



5 < V tracks < 20 
0.1 < Pvis/^s < 0.3 
|pi/ptot| < 0.85 
22 < pt < 50 GeV 
0.55 < T < 0.90 

N icts > 2 
E jet < 25 GeV 
P c > 0.6 



> 10000 GeV 



signal efficiency 



0.340 



0.212 



50 



Table 4: Numbers of generated 
events, and expected events for the 
sequential-cut analysis at \fs = 260 
and 500 GeV for total luminosities of 
50 fb" 1 and 200 ftT 1 with unpolar- 
ized and polarized beams. 







260 GeV 




500 GeV 






generated 


£ = 


501b- 1 


generated 


C = 200 ft)" 1 


P(e-)/P(e + ) 




0/0 


.8/-.6 




0/0 


.8/-.6 


tit* 


50,000 


544 


1309 


50,000 


5170 
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12.4 Iterative discriminant analysis 

The Iterative Discriminant Analysis (IDA) 
[149] is applied to increase the discrim- 
inant power between signal and back- 
ground compared to the sequential-cut- 
based analysis, and thus reduce the sta- 
tistical uncertainty in the stop mass mea- 
surement. Figure 2 gives the results of ex- 
pected number of background events as 
a function of the signal efficiency. The 
chosen working points have efficiencies 
of 38.7% and 41.6% for the y/s = 260 
and 500 GeV analyses, respectively. Ta- 
ble 5 lists the corresponding expected 
background events. 
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Table 5: Numbers of generated events, and expected events for the IDA at 
= 260 and 500 GeV for total luminosities of 50 fb -1 and 200 fb -1 with 
unpolarized and polarized beams. 




Efficiency Efficiency 



Fig. 2: IDA: Expected background events as a function of the signal efficiency. Left: C = 50 fb 1 at -y/s = 260 GeV. Right: 
C = 200 fb" 1 at y/s = 500 GeV. 
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12.5 Systematic uncertainties 

Both the sequential-cut-based analysis and the IDA method lead to a small statistical uncertainty resulting in 
Amjj < 0.2 GeV and thus systematic uncertainties are particularly important to evaluate. Three classes of sys- 
tematic uncertainties are distinguished: 

• instrumental uncertainties related to the detector and accelerator: detector calibration (energy scale), track 
reconstruction efficiency, charm-quark tagging efficiency, and integrated luminosity. 

• Monte Carlo modeling uncertainty of the signal: charm and stop fragmentation effects. The Peterson frag- 
mentation function [150] was used with e c = -0.031 ± 0.011 (OPAL) [151]. For e b = -0.0041 ± 0.0004 
(OPAL) [152] and e b = -0.0031 ± 0.0006 (ALEPH) [153] an average uncertainty of 15% was taken, 
and a factor 2 improvement at the ILC has been assumed, leading to Ae^ = 0.6 x 10~ 6 where ei = 
€b(mh/m mi ) 2 [150, 154]. Fragmentation effects and gluon radiation increase the number of jets signifi- 
cantly and the importance of c-quark tagging is stressed in order to resolve the combinatorics. 

• neutralino mass 108.2 ± 0.3 GeV [155]. 

• theoretical uncertainties on the signal and background. Some improvement compared to the current loop 
calculation techniques is assumed, and an even larger reduction of this uncertainty is anticipated before the 
start of the ILC operation. 

Tables 6 and 7 list the systematic uncertainties for the sequential-cut analysis and the IDA. The systematic 
uncertainty using the IDA method from detector calibration (energy scale) is larger. This is because the sequential- 
cut analysis pays particular attention to cancellation of this uncertainty between the two analyses at the different 
center-of-mass energies. 



Table 6: Sequential-cut analysis experimental sys- 
tematic uncertainties on the signal efficiency. The 
first column indicates the variable that is cut on. 
The second column contains the expected system- 
atic uncertainty for this variable based on experi- 
ence from LEP. The third column shows by how 
much the signal efficiency for y/s = 260 GeV 
varies as a result of varying the cut value by this 
uncertainty. The fourth column shows the same for 
y/s = 500 GeV, and the fifth column lists the re- 
sulting error on the observable Y. 



variable 



error on rel. shift on signal eff. 

variable 260 GeV 500 GeV error on Y 



energy scale 


1% 


3.7% 


3.1% 


0.6% 


N tracks 


0.5% 




negligible 




charm tagging 


0.5% 




taken to be 0.5% 




luminosity 




0.4% 


0.2% 


0.4% 


charm fragmentation 


0.011 


0.3% 


0.8% 


0.6% 


stop fragmentation 


0.6 x 10" 


" 6 0.6% 


0.2% 


0.7% 


neutralino mass 


0.3 GeV 


3.8% 


3.0% 


0.8% 


background estimate 




0.8% 


0.1% 


0.8% 



Table 7: IDA experimental systematic uncertain- 
ties. 
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0.4% 
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0.4% 


charm fragmentation 
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0.1% 


0.6% 


0.5% 


stop fragmentation 


0.6 x 10~' 
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0.8% 


0.7% 


neutralino mass 


0.3 GeV 


3.7% 


1.6% 


2.2% 


background estimate 




0.3% 


0.2% 


0.1% 
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12.6 Mass determination 

The assessment of the achievable stop mass precision is based on the statistical and systematic uncertainties on the 
observable Y (eq. (1)) as summarized in Table 8. The IDA method has a smaller statistical uncertainty, and also a 
smaller background uncertainty due to a smaller number of expected background events. The expected stop mass 
uncertainty is inferred from the uncertainty on Y as given in Table 9. 



Table 8: Summary of statistical and systematic un- 
certainties on the observable Y. 



error source for Y 


sequential cuts 


IDA method 


statistical 


3.1% 


2.7% 


detector effects 


0.9% 


2.4% 


charm fragmentation 


0.6% 


0.5% 


stop fragmentation 


0.7% 


0.7% 


neutralino mass 


0.8% 


2.2% 


background estimate 


0.8% 


0.1% 


sum of experimental systematics 


1.7% 


3.4% 


sum of experimental errors 


3.5% 


4.3% 


theory for signal cross-section 


5.5% 


5.5% 


total error AY 


6.5% 


7.0% 



measurement error Am^ (GeV) 



error category sequential cuts IDA method 
statistical 539 017 
Table 9: Estimated measurement errors (in GeV) sum of experimental systematics 0.10 0.21 
on the stop quark mass. beam spectrum and calibration 0.1 0.1 

sum of experimental errors 0.24 0.28 

sum of all exp. and th. errors 0.42 0.44 



12.7 Cold dark matter interpretation 

The chosen benchmark parameters are compatible with the mechanism of electroweak baryogenesis [141]. They 
correspond to a value for the dark matter relic abundance within the WMAP bounds, ficDM^ 2 = 0.109 [36]. The 
relic dark matter density has been computed as in Ref. [141] 2 . In the investigated scenario, the stop and lightest 
neutralino masses are = 122.5 GeV and m^o = 107.2 GeV, and the stop mixing angle is cos6>j = 0.0105, 
i.e. the light stop is almost completely right-chiral. The improvement compared to Ref. [141] regarding the CDM 
precision determination is shown in Fig. 3 and summarized in Table 10. 



0.14 



0.12 



0.1 



0.08 




121 



122 



123 



124 



Fig. 3: Computation of dark matter relic abun- 
dance QcdmIi 2 taking into account estimated ex- 
perimental errors for stop, chargino, neutralino 
sector measurements at the future ILC. The black 
dots correspond to a scan over the la (Ax 2 < 1) 
region including the total expected experimental 
uncertainties (detector and simulation), the grey- 
dotted region includes also the theory uncertainty, 
and the light grey-dotted area are the previous 
results [141]. The red star indicates the best-fit 
point. The horizontal shaded bands show the la 
and 2a constraints on the relic density measured 
by WMAP. 



2 The assumed benchmark parameters changed slighty (larger slepton masses assumed) and JIcdm/i 2 changed from 
0.1122 [141] to 0.109. 
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Am- tl (GeV) (] C dm/i 2 



Table 10: Estimated precision for the determination of 
stop mass and dark matter relic density for different as- 
sumptions about the systematic errors. 



exp. and th. errors 
stat. and exp. errors only 
sequential-cut analysis 



0.42 o.io9i;;;' 



0.24 



0.109 



.+0.012 
-0.010 
.+0.012 
-0.010 



IDA 



0.28 



0.109 



Conclusions 



Scalar top quarks could be studied with precision at a future International Linear Collider (ILC). The simulations 
for small stop-neutralino mass difference are motivated by cosmology. The precision mass determination at the 
future ILC is possible with a method using two center-of-mass energies, e.g. = 260 and 500 GeV. This method 
can also be applied to other analyses to improve the mass resolution in searches for new particles. The precision 
of two independent analysis methods, one with a sequential-cuts and the other with an Interative Discriminant 
Analysis (IDA) lead to very similar results. The new proposed method increases the mass precision by about 
a factor of three due to the error cancellation using two center-of-mass energies with one near the production 
threshold. Including experimental and theoretical uncertainties, the mass of a 122.5 GeV scalar top could be 
determined with a precision of 0.42 GeV. The interpretation of this benchmark scenario leads to a uncertainty on 
Ocdm'* 2 of —0.013 and +0.015 which is about a factor two better compared to previous results, and comparable 
to current cosmological (WMAP) measurement uncertainties. With the new stop mass determination, the stop 
mass uncertainty is no longer the dominant uncertainty in the QcdmIi 2 calculation. 
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13. Comparison of SUSY spectrum codes: the NUHM case 1 

13.1 Introduction 

Recent analyses of uncertainties in SUSY spectrum calculations [156, 157] have triggered important improvements 
in the various spectrum codes. This concerns in particular the treatment of the top and bottom Yukawa couplings, 
and the finite corrections in translating DR parameters to on-shell masses and mixing angles, see e.g. [158]. 
Moreover, all public codes now apply full two-loop renormalization group (RG) evolution of the SUSY-breaking 
parameters, plus one-loop self-energy corrections for sparticle masses. 

So far, comparisons of spectrum computations have concentrated on the constrained MSSM or mSUGRA 
models. In this contribution, we extend these studies to models with non-universal Higgs masses (NUHM). We 
use the most recent versions of the four public spectrum codes, ISAJET7.75 [159, 158], SOFTSUSY2.0.14 [97], 
SPHEN02.2.3 [18] and SUSPECT 2.3 [37]. We first compare the results from the four public codes for the 
NUHM benchmark points of [160], hence also comparing with the private code SSARD. It turns out that there is 
good agreement once the different sign convention in SSARD as compared to the other codes is taken into account. 
Then we discuss the case of gaugino mediation, for which 10% differences in left-chiral slepton masses can occur 
for very large m 2 Hd — m 2 Hu . 

13.2 NUHM benchmarks 

NUHM models have recently become very popular because they lead to very interesting phenomenological effects 
beyond the well-studied CMSSM case, see e.g. [161, 162, 163]. In general, the non-universality of the Higgs 
doublets can be specified either through GUT-scale masses m 2 ^ and m% d , or though \i and vtia at the weak 
scale. 2 In particular, the group around Baer uses m 2 H , while the group around Ellis and Olive follows the second 
approach. From the four public spectrum codes, only two currently have the (/i, tua) approach implemented. We 
therefore start our discussion by reproducing the NUHM benchmark points proposed in [160] with a scan in the 

TO J7 U _ TO i?d Pl ane - 

Figure 1 shows contours of constant /i and tua from the four public spectrum codes in the plane of m,H u 
versus iriH d , mH u d = sign(m^ u d )\J\ m2 H J' f° r tne NUHM benchmark scenarios a, (5, 7 of [160]. Each 
benchmark point is approximately reproduced where the respective contours of \i and tua intersect. As can be 
seen, the solutions for the various codes lie close to each other. One can therefore also expect good agreement on 
the resulting physical spectrum. 

That this is indeed the case is examplified in Table 1 for benchmark point a. We find agreement of 
about 1-2% on the sparticle and Higgs masses, and about 10% on the relic density, which we compute using 
MICROMEGAS [98, 20] for the ISAJET, SOFTSUSY, SPHENO and SUSPECT spectra. Note also that the agree- 
ment within the public codes themselves is better than the overall agreement including SSARD, cf. Table 1. For 
point (3, there is also 1-2% agreement on the masses. For the relic density, however, the spectra of the public codes 
interfaced to MICROMEGAS give Vth 2 ~ 0.07 with 15% variation, while SSARD gives flh 2 = 0.1. The source 
of the difference is the A-funnel annihilation cross section. 

In SOFTSUSY, there is also the option to input \i and vtia instead of m% at Mqut- To this aim, the 
progam makes an initial guess of the GUT-scale m 2 Hu d for the first iteration. For later iterations, only the EWSB- 
scale boundary conditions are used for m 2 Hu . The procedure works extremely well, giving exactly the same 
results irrespective of whether one uses GUT-scale m 2 Hu or EWSB-scale fi, tua as inputs — at least for the cases 
we have tried. 

There is an analogous option in ISAJET, through the NUSUG keywords. However, the procedure applied 
is more complicated, using boundary conditions for m 2 H at both the GUT and the EWSB scales. The results 
slightly depend on what kind of input is used. For example, for point a, input of \i = 375 GeV and tua — 265 GeV 
gives m 2 Hu = (250) 2 GeV 2 and m 2 Hd = -(320.2) 2 GeV 2 at M GUT . On the other hand, GUT-scale input of 
m 2 Hu = (250) 2 GeV 2 and m 2 Hd = -(320.2) 2 GeV 2 gives (j, = 378 GeV and m A = 274.6 GeV at the weak scale. 

SPHENO and SUSPECT do not yet have n and m A input for SUGRA scenarios. 

'S. Kraml and S. Sekmen 

2 We take pi = /i(Mewsb), where Mewsb is the scale of electroweak symmetry breaking, and tua = mA(pole). 
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NUHM pointy 




mH u 

Fig. 1: Contours of constant /i (solid lines) and uia (dashed lines) from the four public spectrum codes in the mu u versus mn d 
plane for the NUHM benchmark scenarios a, f3, 7 of [160]. The black, red, green and blue lines are for ISAJET, SOFTSUSY, 
SPHENO and SUSPECT, respectively. Each benchmark point is approximately reproduced where the respective contours of p 
and uia intersect. In the grey areas no solution is obtained. 
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SSARD [160] 
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643.9 
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1.79 
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1.46 
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*1 
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1.28 


*2 
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0.50 


h 


590 


599.7 
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594.4 


597.8 


1.63 
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637.4 


630.5 


637.0 


635.4 


1.33 




216 


218.8 


219.3 


218.9 


218.1 


1.51 




296 


296.9 


296.7 


296.7 


295.1 


0.61 




285 


285.1 


285.2 


285.7 


284.6 


0.39 
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298.3 
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0.50 


V T 


285 


284.0 


284.9 


285.4 


284.3 


0.49 


x\ 


115 


112.4 


111.7 


111.9 


112.2 


1.16 




212 


208.2 


207.8 


208.1 


208.2 


2.01 


xi 


388 


380.3 


383.8 


381.7 


380.6 


2.01 


xl 


406 


401.3 


401.5 


400.7 


400.9 


1.32 


xt 


212 


208.3 


209.0 


207.7 


207.6 


2.11 


xt 


408 


400.4 


398.1 


401.5 


400.9 


2.45 


9 


674 


691.2 


687.7 


685.6 


688.7 


2.51 


nti 2 


0.12 


0.109 


0.114 


0.111 


0.108 


10.68 



Table 1: Comparison of results for NUHM benchmark point a given by mo = 210 GeV, rn 1 / 2 = 285 GeV, Aq = 0, 
tan P = 10, m t = 178 GeV, fj, = 375 GeV and m A = 265 GeV. 
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Fig. 2: Neutralino, stau and tau-sneutrino masses (in GeV) from the four public spectrum codes as a function of rriH d , for 
mi/2 = 500 GeV, tan/3 = 10, mO = Ao = m Hu = 500. The black, red, green and blue lines are for ISAJET, SOFTSUSY, 
SPHENO and SUSPECT, respectively. 



13.3 Slepton masses in gaugino mediation 

In general, in SUSY-breaking models with universal scalar and gaugino masses, the right-chiral charged sleptons, 
£r, are lighter than the left-chiral ones and the sneutrinos, £l and i>i (£ = e, fx, r). However, owing to 5-term 
contributions in the RG evolution, for large enough rr? H — rr? H > 0, £l and/or vi can become lighter than £r, 
and even lighter than the Xi [161, 162, 163]. In such a setup, if R parity is conserved, the lightest SUSY particle 
(LSP) should be a gravitino or axino, and the next-to-lightest one (NLSP) a f\ or v r . This has recently attracted 
quite some interest in the context of gaugino mediation [164, 165, 166, 167]. 

In gaugino mediation [168, 169], gauge couplings and gaugino masses each unify at the compactification 
scale Mc, while there are no-scale boundary conditions for sfermion masses and trilinear couplings, i.e. mo = 
A) = 0) [169]. The free parameters of the model are hence m 1( / 2 , m? H , m 2 Hd , tan/3, and the sign of /x; being 
determined by radiative electroweak symmetry breaking. Following [164, 165, 166, 167], we take Mc — Mqut 
and m t = 172.5 GeV. 

Figure 2 shows the Xi< ^i 2 an d masses as a function of mji d , for mw2 = 500 GeV, tan/3 = 10, 
rtiH u = 500. 3 As can be seen, for mn d ~ 1.2 TeV, the left-ciral sleptons become lighter than the right-chiral ones; 
for m,H d ~ 1.6 TeV, they are lighter than the lightest neutralino. Moreover, there is an overall good agreement 
between the codes. Only for very large m? H — m 2 H , when there are large S-term corrections, the differences in 
the left-chiral slepton masses reach ~ 10%. This is also the region where the stau or sneutrino is the (pseudo)LSP. 
Explicit numbers are given in Tables 2 and 3 for m# d = 900 GeV and 1.8 TeV, respectively. Note also that there is 
~ 3-4% difference in fi. The masses of coloured sparticles and Higgs bosons are not shown due to lack of space; 
they agree to ~ 1% or better. 



13.4 Conclusions 

We have compared results of the latest versions of the public SUSY spectrum codes ISAJET, SOFTSUSY, SPHENO 
and SUSPECT for models with non-universal Higgs masses. We find very good agreement for the mass spectra 
and the resulting Vlti 1 . Only for edges of parameter space, the differences in the most sensitive masses become 
large. For instance, we have found differences in slepton masses of the order of 10% for large rr? H — rrijj . 

3 In SUSPECT we use mo = 10~ 3 , since this code does not give a spectrum for mo = 0. 
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2.34 
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301.5 


1.61 
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241.1 


240.9 
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1.29 


Tl 


313.7 


315.7 


315.2 


310.0 


1.82 


V T 


296.6 


300.4 


301.6 


296.7 


1.67 


n — 

x\ 


203.0 


201.8 


202.6 


202.8 


0.59 


xl 


364.4 


365.6 


367.8 


367.2 


0.93 


xl 


463.5 


477.3 


481.2 


477.7 


3.73 


xl 


500.2 


508.4 


513.4 


511.1 


2.6 


xf 


365.0 


367.3 


367.3 


366.3 


0.63 


xf 


499.1 


504.9 


513.1 


510.2 


2.76 



Table 2: Comparison of neutralino, chargino, and slepton masses for a gaugino-mediation scenario with 
mi/j = 500 GeV, tan = 10, m Hu = 500 and m Hd = 900 GeV. 





ISAJET 


SOFTSUSY 


SPHENO 


SUSPECT 


Am [%] 


eR 


432.1 


436.2 


433.8 


436.0 


0.94 


eL 


199.4 


182.7 


192.3 


180.1 


10.23 


Ve 


170.3 


179.9 


175.0 


162.5 


10.12 


Tl 


159.9 


151.1 


163.4 


148.6 


9.5 


T2 


410.5 


413.3 


411.1 


413.0 


0.68 


V T 


130.6 


149.5 


144.5 


148.6 


13.19 


x\ 


202.5 


200.9 


201.3 


200.1 


1.19 




351.7 


353.6 


357.2 


356.4 


1.55 


x° 3 


425.3 


439.4 


444.5 


440.9 


4.39 


xl 


474.9 


481.1 


486.7 


484.6 


2.45 


xf 


351.8 


355.6 


356.2 


355.0 


1.24 


xf 


473.0 


478.9 


486.3 


483.6 


2.77 



Table 3: Same as Table 2 but for m Hd = 1.8 TeV. 
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14. The SLHA2 Conventions 1 
14.1 Introduction 

The states and couplings appearing in the general minimal supersymmetric standard model (MSSM) can be defined 
in a number of ways. In principle, this is not a problem, as translations between different conventions can usually 
be carried out without ambiguity. From the point of view of practical application, however, such translations are, 
at best, tedious, and at worst they introduce an unnecessary possibility for error. To deal with this problem, and to 
create a more transparent situation for non-experts, the original SUSY Les Houches Accord (SLHA1) was proposed 
[72]. However, SLHA1 was designed exclusively with the MSSM with real parameters and R-parity conservation 
in mind. We here summarize conventions [170] relevant for R-parity violation (RPV), flavour violation, and CP- 
violation (CPV) in the minimal supersymmetric standard model. We also consider next-to-minimal models which 
we shall collectively label NMSSM. For simplicity, we still limit the scope of the SLHA2 in two regards: for the 
MSSM, we restrict our attention to either CPV or RPV, but not both. For the NMSSM, we define one catch-all 
model and extend the SLHA1 mixing only to include the new states, with CP, R-parity, and flavour still assumed 
conserved. For brevity, this document only describes our convention choices and not the full ASCII data structures 
that go with them (these are the focus of a complementary summary [171]). The complete SLHA2 is described in 
detail in [170]. 



14.2 Extensions of SLHA1 

Firstly, we allow for using either the A or H + pole masses, respectively, as input instead of the parameter 
m\ (Mi n p U t) defined in [72]. Secondly, we also optionally allow for different parameters to be defined at dif- 
ferent scales (e.g., /i defined at Mewsb, the remaining parameters defined at Mi nput ). 

To define the general properties of the model, we introduce new global switches for field content (either 
MSSM or NMSSM), RPV (either off or on), CPV (either no CPV, just the CKM phase, or general CPV), and 
flavour violation (either no flavour violation or quark and/or lepton flavour violation). 

Also note the recent proposal [172] for a joint SLHA+LHEF (Les Houches Event File [173]) format for 
BSM event generation. 



14.3 Flavour violation 

In the Super-CKM basis of the squarks [174], the quark mass matrix is diagonal, and the squarks are rotated in 
parallel to their superpartners. This implies that only physically measurable parameters are present. Actually, once 
the electroweak symmetry is broken, a rotation in flavour space 

D° = V d D, U° = V U U, D°=U* d D, U°=U* U U, (1) 

of all matter superfields in the (s)quark superpotential 



W Q = e ab (Y D ) ii H?Q b l °D°+(Y u ) ll H b 2 QrU° 



(2) 



brings fermions from the interaction eigenstate basis {d° L , u° L ,d° R , u° R } to their mass eigenstate basis {d^, ul, c£r, ur}: 
d° L = Vddh, u° L = V u ul, d° R — Uddn, u° R = U u ur, and the scalar superpartners to the basis {dL, ul, da, ur}. 
Through this rotation, the Yukawa matrices Yd and Yjj are reduced to their diagonal form Yd and Yjj\ 

{Y D ) lt = {U\YEV d ) u = , {Yu)u = {UlY^V u ) u = . (3) 

V 1 V2 

Tree-level mixing terms among quarks of different generations are due to the misalignment of Vd and V u , expressed 
via the CKM matrix [175, 176], Vckm = V$ Vd, which is proportional to the tree-level ULidLjW + , ULidnjH + , 
and u Ri d Lj H + couplings (i, j = 1, 2, 3). 

In the super-CKM basis the 6 x 6 mass matrices for the up- and down-type squarks are defined as 

£ mass = -$t J M?$ u -$3iM|$ <I , (4) 



'participating/corresponding authors: B.C. Allanach, P. Skands and P. Slavich 
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where = (u L , c L , t L , u R , c R , t R ) T and $ d = (dL,S Ll b L ,d R ,s R ,b R ) T . We diagonalise the squark mass 
matrices via 6 x 6 unitary matrices R u ^d, such that Ru^M'l ji?^ d are diagonal matrices with increasing mass 
squared values. The flavour-mixed mass matrices read: 

Ml ( Vckm roj y c f KM +ml + D u LL - fim u cot \ 

\ ^Tb- - /z*m u cot /3 rh\ + ml + D uRR J 

d V ^T D - fj,*m d t<m(3 m} + m 2 d + D dRR )' 

In the equations above we introduced the 3x3 matrices 

= <J=<X/X«i. ^EC/tTjK, t D = U\TEV A , (7) 

where the un-hatted mass matrices tuq d and trilinear interaction matrices Tjj^d are given in the electroweak 
basis of [72], in which the soft SUSY-breaking potentials V 3 and V 2 have the following forms: 



V 3 = e ab £ [(T B ) y fff + (T D ) l3 H?Qld; R + {Tu) l3 H b 2 Q a lL u* R 



+ h.c. , (8) 



m 4h (to?) 13 u* r + J ifi (m|)y J* R + e iR (mg)ye* R - (mle ab H^H^ + h.c.) . (9) 



The matrices m U; d are the diagonal up-type and down-type quark masses and Df ll,rr are the D-terms 
given by: 

D f ll,rr = cos 2/3 ra| (T| - Q f sin 2 0w) 11 3 , (10) 

which are also flavour diagonal, and Qf is the electric charge of the left-handed chiral supermultiplet to which the 
squark belongs, i.e., it is 2/3 for U and —2/3 for U c . Note that the up-type and down-type squark mass matrices 

in eqs. (5) and (6) cannot be simultaneously flavour-diagonal unless m 2 - is flavour-universal. 

Q 

For the lepton sector, we adopt a super-PMNS basis and cover all cases that lead to a low energy effective 
field theory with Majorana neutrino masses and one sneutrino per family. In terms of this low energy effective the- 
ory, the lepton mixing phenomenon is analogous to the quark mixing case. After electroweak symmetry breaking, 
the neutrino sector of the MSSM contains the Lagrangian pieces 

C = ~v aT (m v )v° +h.c, (11) 

where m v is a 3 x 3 symmetric matrix. The interaction eigenstate basis neutrino fields v° are related to the mass 
eigenstate ones v by v° = V v u, reducing the mass matrix m v to its diagonal form rh v 

{rn v )ii = {Vjm v V v )ii = m Vi . (12) 

The charged lepton fields have a 3 x 3 Yukawa coupling matrix defined in the superpotential piece [72] 

W E = e ab (Y E ) l3 H«L>> E°, (13) 

where the charged lepton interaction eigenstates {e° L , e° R } are related to the mass eigenstates by e° L = V e eL and 
e° R = U e e R . The equivalent diagonalised charged lepton Yukawa matrix is 

(Ye)u = (UlYgVeU = V2^ . (14) 

Lepton mixing in the charged current interaction can then be characterised by the PMNS matrix [177, 178], 
Upmns — VjV v , which is proportional to the tree-level eLiVjW~ and & R iVjH~ couplings. 

Rotating the interaction eigenstates of the sleptons identically to their leptonic counterparts, we obtain 
the super-PMNS basis for the charged sleptons and the sneutrinos, described by the Lagrangian (neglecting the 
possible term ^M%^ v ) 

C rnass = _^ M j^ e _ (15) 
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where = (D e , v^, v r ) T and $ e = (e L , P-L,f L ,e R , ji R , f R ) T . M\ is the 6 x 6 matrix 

,2 , ^,2 [ 7~) fl 'f't 

~2 J £ 
,2 i ^,2 



/ m|.+m^ + C eii ~^p2 e ~ f im e tan /3 \ 



^T E - fi*m e tan /? mj + mj + D eRR 
and M% is the 3 x 3 matrix 

Ml = U PMNS m\UpMNS + D vLL , (17) 
where D eLL and Z?„ LL are given in eq. (10). In the equations above we introduced the 3 x 3 matrices 



in - T 



2 f =V e ^m 2 f V e , ml = Ulmfu e , f E = U^ e T^V e , (18) 



where the un-hatted matrices m| e and Te are given in the interaction basis of ref. [72]. We diagonalise the 
charged slepton and sneutrino mass matrices via the unitary 6x6 and 3x3 matrices i? ei „ respectively. Thus, 
Re,vM~ p Rl v are diagonal with increasing entries toward the bottom right of each matrix. 



14.4 R-parity violation 

We write the R-parity violating superpotential in the interaction basis as 



Wrpv — e ab 



-\ ijk L'?L b j E k + \' ijk L?Qf /•>/..,. - KiL^H^ 



+ \K jk e xyz U?D)Dl (19) 



where x,y,z = 1, . . . , 3 are fundamental SU(3)c indices and e xyz is the antisymmetric tensor in 3 dimensions 
with e 123 = +1. In eq. (19), Xij k , A^- fe and break lepton number, whereas \'!- k violate baryon number. To 
ensure proton stability, either lepton number conservation or baryon number conservation is usually still assumed, 
resulting in either \ ljk = X' ijk = m = or X" jk = for all i,j,k=l, 2, 3. 

The trilinear R-parity violating terms in the soft SUSY-breaking potential are 



^3,rpv — e a fc 



-z{T)ij k L1 L L h jL e* kR + {T')i kL1 L Q b j L d* kR 



+ \(T'% k e xyz ur R q R di R + h.c. . (20) 



Note that we do not factor out the A couplings (e.g. as in T^ k /Xij k = A\^j k ) in order to avoid potential problems 
with \ij k = but Tijk ^ 0. This usage is consistent with the convention for the R-conserving sector elsewhere in 
this report. 

The bilinear R-parity violating soft terms (all lepton number violating) are 

*W = -eabD^Hl + l\ aL m\ iH Hl + h.c. . (21) 

When lepton number is not conserved the sneutrinos may acquire vacuum expectation values (VEVs) 
(^e,Ai,r) = v e,fi,r/ V%- We generalise the SLHA1 parameter v to: 



v\ + v\ + v 2 e +vl + vl = 246GeV. (22) 

The addition of sneutrino VEVs allows for various different definitions of tan (3, but we here choose to keep the 
SLHA1 definition tan (3 — v 2 /vi. 

We use the super-CKM/PMNS basis throughout, as defined in subsection 14.3, with the following consid- 
erations specific to the R-parity violating case. Firstly, the rf-quark mass matrices are given by 

= (YD)i ] v 1 + X kij v k ■ (23) 

where v k are the sneutrino VEVs. Secondly, in the lepton number violating case, we restrict our attention to 
scenarios in which there are no right-handed neutrinos and, thus, neutrino masses are generated solely by the RPV 
couplings. In this case, the PMNS matrix is not an independent input but an output. 

We define the super-CKM basis as the one where the Yukawa couplings Y D and Y v are diagonal. The 
PMNS basis is defined as the basis where Y E is diagonal and the loop-induced neutrino mass matrix is diagonalised. 
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In this way one obtains a uniquely denned set of parameters: 

^ijk = ^rstVv,riVe,sjUe , (24) 
Kjk = Kst V v,riVd,sjUl tk , (25) 

h{ = K r V e r i , (26) 
K 3 k = KstU u , ri U d sj U d j k , (27) 

where the fermion mixing matrices are defined in subsection 14.3. The Lagrangian for the quark-slepton interac- 
tions then takes the following form: 

C = -\' ijk Vid Rk d Lj + Ksk U PMNS,ri V CKM,sjlL,idRkU L j + h.c. . (28) 

Similarly one obtains the soft SUSY breaking couplings in this basis by replacing the superpotential quantities in 
eqs. (24)-(27) by the corresponding soft SUSY breaking couplings. In addition we define: 

As for the R-conserving MSSM, the bilinear terms (both SUSY-breaking and SUSY-respecting ones, in- 
cluding /i) and the VEVs are not independent parameters. Specifically, out of the 12 parameters Ki, Di, sneutrino 
vevs, and m 2 ~ , only 9 are independent. 

Particle mixing 

In general, the neutrinos mix with the neutralinos. This requires a change in the definition of the 4x4 neutralino 
mixing matrix TV to a 7 x 7 matrix. The Lagrangian contains the (symmetric) neutrino/neutralino mass matrix as 

C^ s = -^ 0T M^° + h.c, (30) 

in the basis of 2-component spinors ip° = (v e , v^, v T , —ib, —iw 3 , hi, h-i) T . We define the unitary 7x7 neu- 
trino/neutralino mixing matrix N, such that: 

-\^ T M^ = ~f T N T N*M^ Njj° , (31) 

where the 7 (2-component) generalised neutrinos x° = (yi, v~j) T are defined strictly mass-ordered. 

In the limit of CP conservation, the default convention is that N be a real matrix and one or more of the 
mass eigenstates may have an apparent negative mass. The minus sign may be removed by phase transformations 
on Xi = v i as explained in SLHA1 [72]. 

Charginos and charged leptons may also mix in the case of L-violation. In a similar spirit to the neutralino 
mixing, we define 

qV ass = -\V T M^+ +h.c, (32) 

in the basis of 2-component spinors V + = (e + , n + , r + , — iw + , /i^) T , = (e~ , ,t~ ,—iw~ ,h~^) T where 
— (w 1 =F w 2 )/V2. Note that, in the limit of no RPV the lepton fields are mass eigenstates. 

We define the unitary 5x5 charged fermion mixing matrices U, V such that: 

-\j>- T M^ + = -^^U^U*M^V^\y^ , (33) 

X T diag(m~ + ) X + 

where the generalised charged leptons x + = ( e i , e 2> e t > e t > e t) ^ defined as strictly mass ordered. In the limit 
of CP conservation, U and V are chosen to be real by default. 

R-parity violation via lepton number violation implies that the sneutrinos can mix with the Higgs bosons. 
In the limit of CP conservation the CP-even (-odd) Higgs bosons mix with real (imaginary) parts of the sneutrinos. 
We write the neutral scalars as (fp = ^/2%t(H®, i?2 > v e , v^, v T ) T , with the mass term 

C = -\^ T M\^\ (34) 
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where Ai'io is a 5 x 5 symmetric mass matrix. We define the orthogonal 5x5 mixing matrix H by 



0Tv>T i 



<2 v>T v> .0 



SOT , 2 . <J>0 

diag(m| ) 

where <f>° = {h\, /i!>, /i-;, /14, are the neutral scalar mass eigenstates in strictly increasing mass orderx. 
We write the neutral pseudo-scalars as <f> = v / 2S(_ff°, H^, v ei v^, & T ) T , with the mass term 



(35) 



C = -^ 0T M%$°, 
where M \ a is a 5 x 5 symmetric mass matrix. We define the 4 x 5 mixing matrix H by 



(36) 



(f> K HA^| ^ W , 



(37) 



where <E> = (A", ^4§, A®) are the pseudoscalar mass eigenstates, again in increasing mass order. 

The charged sleptons and charged Higgs bosons also mix in the 8 x 8 mass squared matrix M^± which we 
diagonalize by a 7 x 8 matrix C: 



C = - (H^,H+,e* Lt ,i* Rj )C[ CMl ± C[C 

dia,g(M 2 ±) 



(38) 



where k,l € {l,2,3}and$+ = $"t = , 

R-parity violation may also generate contributions to down-squark mixing via additional left-right mixing 

terms, 



V2 



(39) 



where are the sneutrino vevs. However, this only mixes the six down-type squarks amongst themselves and so 
is identical to the effects of flavour mixing. This is covered in subsection 14.3. 



14.5 CP violation 

In general, we write complex parameters in terms of their real and imaginary parts, rather than in terms of phase 
and modulus. (The SLHA1 data structures are then understood to refer to the real parts, and the imaginary parts 
are provided in data blocks of the same name but prefaced by IM.) The defaults for all imaginary parameters will 
be zero. 

One special case is the fi parameter. When is determined by the conditions for electroweak symmetry 
breaking, only the phase ip^ is taken as an input parameter, see [170]. 

When CP symmetry is broken, quantum corrections cause mixing between the CP-even and CP-odd Higgs 
states. Writing the neutral scalar interaction eigenstates as cfp = \/2(3?^i , ^H®, ^sH®, ^sH®) 7, we define the 
3x4 mixing matrix S by 

- ( f> OT M 2 4>o 0° = -^ T S T S*MloS^ S<j>° , (40) 
^^^^ 

where $° = (h^h^, h®) T are the mass eigenstates. 

For the neutralino and chargino mixing matrices, the default convention in SLHA1 is that they be real. One 
or more mass eigenvalues may then have an apparent negative sign, which can be removed by a phase transfor- 
mation on Xi as explained in SLHA1 [72]. When going to CPV, the reason for introducing the negative-mass 
convention in the first place, namely maintaining the mixing matrices strictly real, disappears. In the CP violating 
case, we therefore take all masses real and positive, with N, U, and V complex. 
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14.6 The next-to-minimal supersymmetric standard model 

We write the most general CP conserving NMSSM superpotential as: 

Wnmssm = Wmssm - e ab \SH^H b 2 + ^ K S 3 + p'S 2 + H F S , (41) 

where Wmssm is the MSSM superpotential, in the conventions of ref. [72, eq. (3)]. A non-zero A in combination 
with a VEV (S) of the singlet generates a contribution to the effective fi term fi e s = A (S) + fi, where the MSSM 
/i term is normally assumed to be zero in NMSSM constructions, yielding /j, e s = A (S). The remaining terms 
represent a general cubic potential for the singlet; k is dimensionless, p! has dimension of mass, and £p has 
dimension of mass squared. The soft SUSY-breaking terms relevant to the NMSSM are 

Koft - V 2 . MSS M + Vzmssm + ™ 2 \S\ 2 + (-e ab \A x S H<£ H b 2 + ^ K A K S 3 + m' 2 S 2 + £ S S + h.c.) , (42) 

where Vi mssm are the MSSM soft terms defined in eqs. (8) and (9), and we have introduced the notation m'g = 
B'n'. 

At tree level, there are thus 15 parameters (in addition to niz which fixes the sum of the squared Higgs 
VEVs) that are relevant for the Higgs sector of the R-parity and CP-conserving NMSSM: 

tan/3, fj,, m 2 Hl , m 2 H2 , m\, A, k, A x , A k , fj,', m'|, £ F , £ s , A (S) , m% . (43) 

The minimisation of the effective potential imposes 3 conditions on these parameters, such that only 12 of them 
can be considered independent. 

Note that we write the soft parameter m% in the form m§/ (cos (5 sin (3), see [72]. The notation m 2 A that was 
used for that parameter in the SLHA1 is no longer relevant in the NMSSM context, but by keeping the definition in 
terms of m 2 and cos (3 sin (3 unchanged, we maintain an economical and straightforward correspondence between 
the two cases. 



Particle mixing 

In the CP-conserving NMSSM, the CP-even interaction eigenstates are <j)° = y/2$t(H$, H$, S) T . We define the 
orthogonal 3x3 mixing matrix S by 

-cjP T Mlo<P° = -^J^SM 2 oS T S£ , (44) 

where $° = (h^h®, ^3) are the mass eigenstates ordered in mass. 

In the CP-odd sector the interaction eigenstates are 0° = \^2^s (Hi , H$, S) T . We define the 2 x 3 mixing 
matrix P by 

-^ T Ml ^°^~^P^PM 2 P T P^, (45) 

where $° = (A®, A®) are the mass eigenstates ordered in mass. 

The neutralino sector of the NMSSM requires a change in the definition of the 4 x 4 neutralino mixing 
matrix to a 5 x 5 matrix. The Lagrangian contains the (symmetric) neutralino mass matrix as 

£ F S = -l$ OT Mjoj>° + h.c. , (46) 

in the basis of 2-component spinors tp = (—ib, —iw 3 , hi,h 2 , s) T . We define the unitary 5x5 neutralino mixing 
matrix N such that: 

-\^° T M^ a - -\i^N^N*M i0 N^ Ag , (47) 

X 0T diag(m-o) ^° 

where the 5 (2-component) neutralinos \% are defined such that the absolute value of their masses increase with i. 
As in SLHA1, our convention is that A be a real matrix. One or more mass eigenvalues may then have an apparent 
negative sign, which can be removed by a phase transformation on \i- 
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14.7 Summary 

The Supersymmetry Les Houches Accord (SLHA) [72] provides a universal set of conventions for supersymmetry 
analysis problems in high energy physics. Here, we summarise extensions of the conventions of the first SLHA 
to include various generalisations [170]: the minimal supersymmetric standard model with flavour violation, RPV, 
and CPV, as well as the simplest next-to-minimal model. For updates and examples, see 

http : / /home . f nal . gov/^skands/ slha/ 
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